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During  the  second  year  (FY'82)  a  great  deal  of  process  has  been  made 
on  the  white-light  optical  information  processing  and  holography  research 
program.  In  this  period,  we  have  evaluated  the  coherence  requirement, 
source  encoding,  and  signal  sampling  concept  for  the  proposed  white-light 
optical  signal  processor.  We  have  shown  that  the  spatial  coherence 
requirement  is  governed  by  source  distribution,  while  the  temporal 
coherence  requirement  is  controlled  by  spatial  frequency  bandwidth  of  the 
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input  signal  and  the  grating  sampling  frequency.  In  order  to  alleviate 
the  basic  constraints  of  a  white-light  source,  we  have  developed  a  source 
encoding  and  signal  sampling  concept,  so  that  the  information  can  be 
processed  in  complex  amplitude  for  the  entire  spectral  band  of  the  white- 
light  source.  We  have  also  evaluated  an  apparent  transfer  function  for 
the  proposed  white-light  signal  processor.  We  have  shown  that  the  MTF 
is  dependent  upon  the  degree  of  spatial  and  temporal  coherence.  The 
derived  apparent  transfer  function  is  very  general,  which  can  be  applied 
to  any  partially  coherent  optical  processor.  Since  the  proposed  white- 
light  signal  processor  is  very  suitable  for  color  image  processing,  we 
have,  in  this  period,  also  demonstrated  several  color  image  processing 
capabilities.  Among  those  are  broadband  color  image  deblurring  and  color 
image  subtraction.  From  those  results,  we  have  seen  high  quality  de- 
blurred  color  images  and  subtracted  color  images  can  be  performed  by  the 
proposed  white-light  optical  processor.  We  have  also  in  this  period 
evaluated  the  primary  aberrations  and  the  bandwidth  requirements  for 
rainbow  holographic  processes.  The  conditions  for  the  removal  of  the 
five  primary  aberrations  are  also  determined. 
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I.  Introduction 

We  have  in  the  past  year  investigated  several  research  projects  on 
"White-Light  Optical  Information  Processing  and  Holography."  Our  research 
works  have  been  quite  consistent  with  our  AFOSR  sponsored  research  grant 
reporting  in  various  journals  and  conference  publications.  Sample 
copies  of  these  papers  are  included  in  this  annual  report  in  the  sub¬ 
sequent  sections,  for  which  to  provide  a  concise  documentation  of  our 
work.  In  the  f ollowing  sections ,  we  will  provide  a  general  overview 
of  our  research  progress  made  in  the  past  twelve  months  and  we  will 
detail  some  of  those  research  progresses.  A  list  of  publications 
resulting  from  AFOSR' 8  support  is  Included  at  the  end  of  this  report. 

II.  Summary  and  Overview 

Since  the  invention  of  laser  (i.e.,  a  strong  coherent  source)  laser 
has  become  a  fashionable  tool  for  many  scientific  applications  particu¬ 
larly  as  applied  to  coherent  optical  signal  processing.  However  coherent 
optical  signal  processing  systems  are  plagued  with  coherent  noises, 
which  frequently  limit  their  processing  capability.  As  noted  by  the 
late  Gabor,  the  Nobel  prize  winner  in  physics  in  1970  for  his  Invention 
of  holography,  the  coherent  noise  is  the  number  one  enemy  of  the  Modern 
Optical  Signal  Processing  [lj.  Aside  the  coherent  noise,  the  coherent 
sources  are  usually  expensive,  and  the  coherent  processing  environments 
are  very  stringent.  For  example,  heavy  optical  benches  and  dust  free 
environments  are  generally  required. 

Recently,  we  have  looked  at  the  optical  processing  from  a  different 


standpoint.  A  question  arises,  is  it  necessarily  true  that  all  optical 
signal  processing  required  a  coherent  source?  The  answer  to  this  question 
is  that  there  are  many  optical  signal  processings  that  can  be  carried  out 
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by  a  white-light  source  [2].  The  advantages  of  the  proposed  white-light 
signal  processing  technique  are:  1.  It  is  capable  of  suppressing  the 
coherent  noise;  2.  White-light  sources  are  usually  inexpensive:  3.  The 
processing  environments  are  not  critical;  4.  The  white-light  system  is 
relatively  easy  and  economical  to  maintain;  and  5.  The  white-light  pro¬ 
cessor  is  particularly  suitable  for  color  image  processing. 

There  is  however  a  basic  different  approach  toward  a  coherent  and 
a  white-light  processor.  In  coherent  processing,  virtually  no  one  Beems 
to  care  about  the  coherence  requirements,  since  the  laser  provides  a 
strong  coherent  source.  However,  in  white-light  processing,  the  knowledge 
of  the  coherence  requirement  is  usually  needed. 

In  white-light  processing  we  would  approach  the  problem  backward. 
First,  we  should  know  what  is  the  processing  operation  we  wish  to  perform: 
Is  it  a  1-D  or  2-D  processing?  Is  the  signal  filtering  a  point  or  point- 
pair  concept?  What  is  the  spatial  bandwidth  of  the  signal?  etc.  Then 
with  these  knowledge,  we  would  be  able  to  evaluate  the  coherence  require¬ 
ments  at  the  Fourier  and  at  the  input  planes.  From  the  evaluated  results, 
we  would  be  able  to  design  a  signal  sampling  function  and  a  source  encoding 
function  to  obtain  these  requirements.  The  objective  of  using  a  signal 
sampling  function  is  to  achieve  a  high  degree  of  temporal  coherence  in 
Fourier  plane  so  that  the  signal  can  be  processing  in  complex  amplitude, 
for  the  entire  spectral  band  of  a  white-light  source.  And  for  the  source 
encoding  is  to  alleviate  the  constraints  of  an  extended  white-light  source. 

In  the  following  sections,  we  shall  highlight  some  of  the  research 
done  during  the  past  year  effort  on  the  white-light  signal  processing. 
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2.1  Coherence  Requirement  (See  Section  III) 

Our  research  In  the  past  year  period  has  included  the  evaluation 
of  the  coherent  requirement  for  white-light  signal  processing  [3],  The 
mutual  intensity  function  for  a  partially  coherent  light  is  used  to 
develop  an  expression  for  the  output  intensity  distribution  for  a  broad¬ 
band  optical  information  processor.  The  coherence  requirement  for 
smeared  image  deblurring  and  image  subtraction  is  then  determined  using 
the  intensity  distribution. 

We  have  shown  that  the  temporal  and  spatial  coherence  requirements 
for  some  partially  coherent  optical  processing  operations,  namely, 
image  deblurring  and  image  subtraction,  can  be  determined  in  terms  of 
the  output  intensity  distribution.  For  image  deblurring  the  temporal 
coherence  requirement  depends  on  the  ratio  of  the  deblurring  width  to 
the  smeared  length  of  the  blurred  image.  To  obtain  a  higher  degree  of 
deblurring  a  narrower  spectral  width  of  the  light  source  is  required. 

For  example,  if  the  deblurring  ratio  Aw/w  is  0.1,  the  spectral  width, 
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AX,  should  be  <  640  A. 

For  the  spatial  coherence  requirement  the  image  deblurring  depends 
on  both  the  deblurring  ratio  Aw/w  and  the  smeared  length  w.  If  the 
deblurring  ratio  Aw/w  -  1/10  and  w  ■  1  mm,  a  slit  source  <  0.26  mm 
should  be  used.  For  a  smeared  image  deblurring  operation  the  con¬ 
straints  of  the  temporal  and  spatial  coherence  requirements  are  not 
critical,  which  can  be  achieved  in  practice. 

For  image  subtraction,  the  temporal  coherence  requirement  is 
determined  by  the  highest  spatial  frequency  and  the  separation  of  the 
input  object  transparencies.  If  the  separation  and  spatial  frequency 
of  the  input  transparencies  are  high,  a  narrower  spatial  bandwidth  of 
the  light  source  is  required. 
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For  spatial  coherence  requirement,  the  modulation  transfer  function, 
which  determines  the  contrast  of  the  subtracted  image,  depends  on  the 
ratio  of  the  slit  width  to  the  spatial  period  of  the  encoding  mask, 
ie.,  d/0.  If  the  ratio  d/D  is  low,  a  higher  contrast  subtracted  image 
can  be  obtained.  For  example,  with  d/D  ■  0.05,  a  relatively  higher 
MTF  *  0.85  can  be  obtained.  Compared  with  the  image  deblurring  operations, 
the  coherence  requirements  are  more  stringent  for  the  subtraction  process. 

Finally,  we  have  concluded  that,  the  solution  to  the  coherence 
requirement  for  partially  coherent  processing  is  not  restricted  to  the 
application  of  the  deblurring  and  subtraction  operation,  but  may  also 
be  applied  to  any  other  optical  processing  operation. 

2.2  Apparent  Transfer  Function  (See  Section  IV). 

We  have  in  the  past  year  period  evaluated  an  apparent  transfer 
function  for  our  white-light  optical  signal  processing  system  [4], 

The  nonlinear  behavior  of  the  partially  coherent  optical  processor, 
when  considering  either  Intensity  or  amplitude  distribution  input  signals, 
necessitates  the  use  of  the  apparent  transfer  function  to  accurately 
predict  the  system  response.  We  have  derived  the  general  formulas  for 
MTF  in  terms  of  the  theory  of  partially  coherent  light.  These  deriva¬ 
tions  indicate  the  dependence  of  MTF  upon  the  degree  of  spatial  cbher- 
ence  (i.e.,  the  source  size)  as  well  as  the  degree  of  temporal  coherence 
(i.e.,  the  source  spectral  bandwidth).  MTF  has  been  shown  to  be  less 
dependent  upon  the  spatial  coherence  requirement  as  compared  to  its 
relationship  with  the  temporal  coherence  requirement. 

It  has  been  oted  that  the  spatial  bandwidth  of  our  optical  processor 
is  primarily  dependent  upon  the  size  of  the  filter  Ax,  where  the  filter 
is  placed  in  the  Fourier  plane.  The  transfer  systems  bandwidth  may  be 


increased  by  using  a  larger  spatial  filter  Ax.  However,  the  size  of 
the  filter  is  selected  such  that  Ax  -  pQfAX,  which  is  linearly  related 
to  the  spatial  frequency  of  grating,  the  focal  length  of  the  transform 
lens  and  the  spectral  width  of  the  light  source.  A  narrow  spectral 
band  AX  is  necessary  for  most  partially  coherent  optical  information 
processing  operations.  In  order  to  achieve  the  required  AX  for  a 
wide  strip  of  spatial  filter  Ax  in  the  spatial  frequency  plane,  a  dif¬ 
fraction  grating  of  sufficiently  high  frequency  pQ  at  the  input  plane 
is  needed.  For  example,  for  partially  coherent  processing  with  a 
white  light  source,  a  set  of  narrow  spectral  band  filters,  each  with 
a  spectral  bandwidth  AX,  can  be  used  in  the  spatial  frequency  plane. 

Finally,  we  conclude  that  the  apparent  transfer  function  which 
we  have  obtained  is  rather  general  and  may  be  applied  to  any  partially 
coherent  optical  processing  system. 

2.3  Source  Encoding  and  Signal  Sampling  (See  Section  V). 

He  have  in  the  past  year  developed  a  concept  of  source  encoding 
and  signal  sampling  for  our  white-light  signal  processor  [5].  We  have 
shown  that  the  advantage  of  source  encoding  is  to  provide  an  appropriate 
spatial  coherence  function  at  the  input  plane  so  that  the  signal  processing 
can  be  carried  out  by  an  extended  incoherent  source.  The  effect  of  the 
signal  sampling  is  to  achieve  the  temporal  coherence  requirement  at  the 
Fourier  plane  so  that  the  signal  can  be  processed  in  complex  amplitude. 

If  the  filtering  operation  is  two-dimensional,  a  multi-spectral-band  2-D 
filters  should  be  utilized.  If  the  filtering  operation  is  one-dimensional, 
a  fan-shape  filter  can  be  used. 

In  short,  one  should  carry  out  the  processing  requirements  backward 
for  a  white-light  processing.  With  these  processing  requirements  (e.g.. 


operation,  temporal  and  spatial  coherence  requirements),  multi-spectral - 
band  or  broad-band  filter,  signal  sampling  function,  and  source  encoding 
mask  can  be  synthesized.  Thus  the  signal  processing  can  be  carried  out 
in  complex  amplitude  over  the  whole-spectral  band  of  the  white-light 
source. 

2.4  Broadband  Image  Deblurring  (See  Reference  6) 

In  this  period  of  research,  we  have  shown  a  broadband  color  image 
deblurring  technique  utilizing  a  white*-light  source  [6].  This  broad 
spatial  band  deblurring  technique  utilized  a  grating  base  method  to 
obtain  a  dispersed  smeared  image  spectra  in  the  Fourier  plane  so  that 
the  deblurring  operation  can  be  taken  placed  in  complex  amplitude  for 
the  entire  visible  wavelengths.  To  perform  this  complex  amplitude 
deblurring  for  the  entire  spectral  band  of  the  light  source^  we  have 
shown  that  a  fan-type  deblurring  filter  to  compensate  the  scale  varia¬ 
tion  of  the  smeared  signal  spectra  due  to  wavelength  dispersion  can  be 
utilized.  To  alleviate  the  low  transmission  efficiency  of  the  deblurring 
filter,  we  synthesized  the  deblurring  filter  with  the  combination  of  a 
broadband  phase  filter  and  a  fan-shaped  amplitude  filter.  The  broad 
spectral  band  phase  filter  is  synthesized  by  optical  coating  techniques, 
whil^  the  fan- shaped  amplitude  filter  is  obtained  by  a  1-D  coherent 
processing  technique. 

By  comparison  of  the  results  obtained  by  the  broadband  image 
deblurring  with  the  narrow  spectral  band  and  coherent  techniques, 
we  have  seen  that  the  results  obtained  by  the  broadband  deblurring  offer 
a  higher  image  quality.  We  have  also  shown  that  the  broadband  deblurring 
technique  is  very  suitable  for  color  image  deblur ring.  We  have  provided 
several  color  image  dsblurrlng  results  Obtained  by  the  broadband  deblurring 
technique.  From  these  color  debVurred  images  we  have  seen  that  the 
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fidelity  of  the  color  reproduction  is  very  high  and  the  quality  of 
deblurred  image  is  rather  good.  Although  there  is  some  degree  of 
color  blur  due  to  chromatic  aberration  of  the  transform  lenses,  it 
can  be  eliminated  by  utilizing  higher-quality  achromatic  transform 
lenses. 

Further  improvements  of  the  deblurring  can  also  be  obtained  by 
utilizing  a  blazed  grating  to  achieve  a  higher  smeared  spectral 
diffraction  efficiency  so  that  a  wider  spatial  band  deblurring  filter 
can  be  used  to  achieve  a  higher  degree  of  deblurring. 

2.5  Color  Image  Subtraction  (See  Section  VI) 

We  have  also  in  this  period  accomplished  a  research  project  on 
color  image  subtraction  with  extended  incoherent  sources  [7].  We 
have  introduced  a  source  encoding  technique  to  obtain  a  point-pair 
spatial  coherence  function  for  the  subtraction  operation.  The  basic 
advantage  of  source  encoding  is  to  Increase  the  available  light  power 
for  the  image  subtraction  operation,  so  that  the  inherent  difficulty 
of  obtaining  Incoherent  point  sources  can  be  alleviated.  Since  the 
technique  uses  incoherent  sources,  the  annoying  coherent  artifact 
noise  can  be  suppressed.  We  would  see  that  the  concept  of  color  image 
subtraction  may  also  be  extended  to  the  use  of  white-light  source,  for 
which  a  program  is  currently  under  investigation.  In  experimental 
demonstrations,  we  have  shown  that  color  subtracted  Images  can  be  easily 
obtained  by  this  incoherent  subtraction  technique.  Since  virtually  all 
images  of  natural  objects  are  color,  the  technique  may  offer  a  wide 
range  of  practical  applications. 

2.6  Rainbow  Holographic  Aberrations  (See  Section  VII) 

In  this  period,  ve  have  also  theoretically  evaluated  the  primary 
aberrations  and  the  bandwidth  requirements  for  rainbow  holographic 


processes  [8] .  The  results  obtained  for  the  rainbow  holographic  process 
are  rather  general,  for  which  the  conventional  holographic  image 
resolution,  aberrations,  and  bandwidth  requirements,  can  be  derived. 

The  conditions  for  the  elimination  of  the  five  primary  rainbow  holo¬ 
graphic  aberrations  are  also  given.  These  conditions  may  be  useful 
for  the  application  of  obtaining  a  high-quality  rainbow  hologram  image. 
In  terms  of  bandwidth  requirements,  we  have  shown  that  the  bandwidth 
requirement  for  a  rainbow  holographic  construction  is  usually  several 
orders  lower  than  that  of  a  conventional  holographic  process.  There¬ 
fore,  a  lower -resolution  recording  medium  can  generally  be  used  for 
most  of  the  rainbow  holographic  constructions. 

2 . 7  Remarks 

We  have  demonstrated  in  the  past  twelve  months  that  the  proposed 
white-light  optical  signal  processor  is  very  suitable  for  color  image 
processings,  which  is  one  step  closer  to  the  reality  of  the  practical 
optical  image  processing.  As  compared  with  the  conventional  coherent 
optical  signal  processor  and  the  electronic  digital  counterparts, 
the  white-light  signal  processor  offers  the  advantage  of 

1.  Cost  effective  in  operation, 

2.  Simplicity  of  system  maintenance, 

and  3.  Suitability  for  color  signal  processing. 

However,  there  is  a  basic  different  approach  toward  a  coherent 
and  a  white-light  processor.  In  white-light  processing,  the  knowledge 
of  the  coherence  requirement  is  usually  needed.  With  the  knowledge  of 
the  temporal  coherence  requirement  in  the  Fourier  plane  and  the  spatial 
coherence  requirament  at  the  input  plane,  a  signal  sampling  function  and 
a  source  encoding  function  can  be  utilized  to  achieve  these  requirements 
such  that  the  signal  processing  can  be  carried  out  by  an  extended  white- 
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light  source.  In  other  words,  one  should  carry  out  the  processing 
requirements  backward  for  a  white-light  processing,  so  that  the 
signal  can  be  processing  in  complex  amplitude  rather  than  in  intensity. 

In  short,  we  have  shown  that  the  white-light  processor  on  one 
hand  can  process  the  signal  in  complex  amplitude  like  a  coherent 
processor,  on  the  other  hand  it  suppresses  the  coherent  noise  like 
an  incoherent  processor. 

2.8  Future  Research 

Our  future  research  is  expected  to  follow  the  general  direction 
addressed  in  the  proposal.  Attention  will  be  given  in  the  following: 

1.  To  develop  a  real-time  processing  capability  for  the 
proposed  white-light  processing  technique. 

2.  To  carry  out  various  experimental  demonstrations  for  the 
real-time  processing  mode. 

3.  To  develop  a  real-time  color  image  processing  technique. 

4.  To  carry  out  the  real-time  color  image  processing  capa¬ 
bilities  as  applied  to  color  image  transformation,  restora¬ 
tion,  detection,  encoding,  etc. 

5.  To  develop  a  technique  of  synthesizing  a  multi-wavelength 
complex  spatial  filter  that  is  suitable  for  white-light 
signal  detection  system. 

6.  To  carry  out  a  computer  generated  spatial  filter  program  for 
the  broadband  optical  detection  and  recognition  problem. 

7.  Provides  experimental  demonstrations  of  the  principles  and 
processing  operation  that  we  propose. 

8.  Extension  of  the  source  encoding  technique  for  image  subtrac¬ 
tion  utilizing  a  broadband  white-light  source. 


ir*i 


9.  The  implementation  of  natural  sun  light  for  the  proposed 
processing  system. 

In  summary,  our  goal  is  to  develop  a  full  research  program  on  the 
white-light  signal  processing  technique  which  has  the  capability  of 
carrying  out  all  the  major  processing  operations  that  a  coherent 
optical  processor  can  offer.  This  white«-light  optical  processing  system 
is  vitally  essential  in  many  areas  of  monochrome  and  color  signal 
processings.  We  are  confident  that  this  white-light  processing 
technique  will  provide  a  step  toward  the  research  and  development 
of  white-light  and  sunlight  optical  computers. 
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Coherence  requirements  for  partially  coherent  optical 
processing 

S.  L.  Zhuang  and  F.  T.  S.  Yu 


The  mutual  intensity  function  for  a  partially  coherent  light  is  used  to  develop  an  expression  for  the  output 
intensity  distribution  for  a  broadband  optical  information  processor.  The  coherence  requirement  for 
smeared  image  deblurring  and  image  subtraction  is  then  determined  using  the  intensity  distribution.  We 
also  quantitatively  show  the  dependence  of  coherence  criteria  on  the  spectral  bandwidth,  the  source  size, 
deblurring  width,  spatial  frequency,  and  the  separation  of  input  object  transparencies. 


I.  Introduction 

Optical  systems  perform  myriad  sophisticated  in¬ 
formation  processing  operations  using  coherent  light.1-3 
However,  coherent  systems  are  susceptible  to  coherent 
artifact  noise,  which  frequently  limits  their  processing 
capability.  The  use  of  incoherent  light  alleviates  the 
noise  problem,4-7  usually  at  the  cost  of  a  lowered  sig- 
nal-to-noise  ratio  due  to  dc  bias.  Various  workers  have 
studied  optical  systems  which  employ  (a)  a  totally  in¬ 
coherent  source8-9  or  (b)  a  broadband  source  with  re¬ 
duced  coherence. 10  The  question  to  be  addressed  in  the 
paper  is,  to  what  degree  can  we  relax  the  coherence  re¬ 
quirement  without  sacrificing  the  overall  results  of  the 
processing  system?  We  use  Wolfs11  theory  of  partially 
coherent  light  to  develop  the  necessary  coherence  cri¬ 
teria  for  an  optical  processor.  The  results  are  used  to 
discuss  the  temporal  and  spatial  coherence  require¬ 
ments  for  the  specific  problems  of  image  deblurring  and 
image  subtraction. 

The  nature  of  optical  processing  operations  governs 
the  spatial  and  temporal  coherence  requirements  nec¬ 
essary  to  obtain  satisfactory  results.  As  system  con¬ 
figuration  varies  from  one  application  to  the  other  it  is 
desirable  to  develop  a  generalized  function  for  the  sys¬ 
tem  output  intensity  expressed  as  a  function  of  the 
spectral  and  spatial  bandwidths. 

Consider  the  partially  coherent  optical  processing 
system12-14  in  Fig.  1.  Let  the  source  plane  Po  be  that 
of  an  incoherent  spatially  extended  light  source  (i.e., 
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white-light  source)  with  an  intensity  distribution 
y(xojo).  The  input  plane  Pi  contains  the  signal 
transparency  and  also  a  coding  mask  (e.g.,  a  sinusoidal 
grating)  if  needed.12-14  Let  the  complex  amplitude 
transmittance  at  this  plane  be  t(u,u).  Its  spectrum  may 
then  be  processed  in  complex  amplitude  by  placing  a 
complex  spatial  filter  f(xy)  in  the  Fourier  plane.  The 
theory  of  partially  coherent  light  can  be  used  now  to 
write  the  output  intensity  distribution,11-16 

XAo+AA/2  - 

-AA/2  JJ  ■),<*o«yo)S<A>C<A> 

*IJT.  T(x  +  xc.v  +  vo;A)/(x,y> 

X  expj-;  ^  (xu'  +  yi/)Jdxdy|2  dx^dygd A,  ( 1 1 

where  S(A)  and  C(X)  are  defined  to  be  the  relative 
spectral  intensity  of  the  light  source  and  the  relative 
spectral  response  sensitivity  of  the  detector,  respec¬ 
tively.  r(x,y;A)  is  the  Fourier  spectrum  of  the  input 
transparency  for  wavelength  A  of  the  light  source.  The 
integration  is  performed  over  the  spectral  bandwidth 
AA  of  the  light  source,  Ao  being  its  center  wavelength. 

In  the  following  sections  we  shall  use  the  output  in¬ 
tensity  distribution  given  by  Eq.  (1)  to  evaluate  the 
coherence  requirement  for  image  deblurring  and  image 
subtraction  of  a  partially  coherent  optical  information 
processor. 

II.  Coherence  Requirement  for  Image  Demurring 

A  photographic  image  debiurring  technique 
employing  a  white-light  source  has  been  described  in 
previous  papers.13-15  We  had  briefly  stated  that  the 
coherence  requirements  depend  on  the  smeared  length 
of  the  object  and  the  size  and  spectral  bandwidth  of  the 
light  source. 17  We  shall  now  discuss  the  temporal  and 
spatial  coherence  requirements  for  image  debiurring 
separately. 
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A.  Temporal  Coherence  Requirement 
We  consider  a  spectrally  broadband  point  source  for 
which  the  intensity  distribution  can  be  expressed  by  a 
Dirac  b  function,  i.e.,  "yUonyo)  ■  5(jto^0)-  The  spectral 
distribution  is  uniform  over  the  bandwidth  [s(X)  *  k, 
a  constant].  If  the  smeared  length  W  is  known  a  priori, 
the  output  deblurred  image  irradiance  can  be  derived 
from  Eq.  (1)  as 


/<u>')-  f  **  |A(u>';X)|*dX,  (2) 

where  X/  and  X*  are  the  extreme  lower  and  upper 
wavelengths  of  the  light  source; 


A(u' »  J*  T(xwv;X)/(x,y)  exp 

is  the  output  complex  light  distribution  of  the  deblurred 
image  due  to  wavelength  X;  T(x,y;\)  is  the  Fourier 
spectrum  of  the  input  blurred  transparency  t(u,v)  due 
to  X;  and  /(x,y)  is  the  deblurring  filter. 

We  note  that,  if  the  input  blurred  image  is  superim¬ 
posed  with  a  sinusoidal  phase  grating,  the  blurred  image 
spectrum  will  disperse  into  rainbow  colors  in  the  Fourier 
plane  thus  allowing  a  stripwise  design  of  deblurring 
niter  for  each  narrow  spatial  band  in  the  Fourier  plane. 
It  is  evident  that  the  temporal  coherence  requirement 
is  limited  by  the  narrow  spectral  band  (Le.,  narrow  strip) 
of  the  deblurring  filter.  If  the  input  object  is  a  linear 
smeared  image,  a  fan-shaped  deblurring  niter  can  be 
used  to  compensate  the  wavelength  variation.13-1'’ 
Since  each  of  the  narrow  strip  deblurring  niters  is  spa¬ 
tial  frequency-limited  over  a  narrow  spectral  band,  the 
coherence  requirement  for  each  narrow  strip  filter  can 
be  equivalently  analyzed  as  the  case  without  the  input 
phase  grating  (i.e.,  the  coding  mask). 

Since  the  image  deblurring  takes  place  at  every  image 
point  the  linear  blurred  image  transparency  may  be 
represented  as 


.  2t 

-j  —ixu  +yv  ) 
X/ 


dxdy  (3) 


Fig.  1.  Partially  coherent  optical  processing  system:  P0,  source 
plane;  Pit  input  plane;  P2,  Fourier  plane;  P3,  output  plane;  L,  achro¬ 
matic  lenses. 


t(u.i->  “  rect(u/W),  (4) 

T(x,y;X)  -  J|r(u,t  )|  “  sine  jj^xj  •  (5) 

where  f  is  the  focal  length  of  the  achromatic  transform 
lenses,  and 


recttn/W)  «  I 
10 

The  deblurring  filter  is  thus 


|u|  «  W/2, 

otherwise. 


(6) 


/(x:X0> 


1 

T(x;X0> 


(7) 


where  Xo  is  the  center  wavelength  of  the  light  source, 
and  for  simplicity  we  have  adopted  a  I-D  notation.  The 
corresponding  deblurred  image  is  given  by 


<4(u';X)  ■  J*  T(x;X)/(x;A„)  exp|-y  ^  xu'Jdx 

i'/W)»  f  - - - expf-y  77xu')dx, 

J—  iw  \  \  x/  ; 


1  recti  u'i 


W_ 
,X„ f 


(81 


where  *  denotes  the  convolution  operation.  By  a 
straightforward  but  tedious  calculation,  illustrated  in 
the  Appendix,  the  deblurred  image  of  Eq.  (8)  can  be 
shown  to  be 


A(u';X)  - 


4X</  iS  ,  ...  .  '"X)  .  2anXo  , 
-  E  <-l)"sin  — —  sin— — -u 

X  n  —  1  X  X IV 


sgnlu'l,  for|uj>W72, 


4X«f  2" 

-  —  I  (-1)" 

X  n  - 1 


xnXo  2  vnKi 
cos — - — cos  ....  u' 


(9) 


X  XW 

where  the  series  converges  for  2m  »  1.  In  the  final 
result  it  can  be  shown  that  the  output  image  irradiance 
is  given  by 

/«>'(u',AX)  E  E  (3|X;,cos|»'.(u')X.l 

r‘  s  s' l-l 

-XI3  cos|al,„  (u')X»||  -  6  (a',.  <u')|  IX;*  sinja!,.  (u')X«| 

-  X;*  sinja !,„  (u')Xsll  -  6  (aL. (u'»|2  |X;‘  cos[a!,„  lu'IX.) 

-  X;1  casjaj,,.  (u')Xsl)  -  8  |a!,„  (u'll3  {SilX,)  -  Sif X* >|>,  forjuj  >  W/2, 


for  |u'|  «  W/2, 


(10) 


/•*'(«',  AX) -^EIZ  <-l>*+"’  <6|o' „  (u')|3  ]Si(X.)  -  Si(X»)| 

X*  s  s'  l-l 

+  6(a!,.  (u')|2|X;’  cosja',,  (u'lX.1  -  X;1  cos|a'„„  <u')X»|| 

+  6loi.  (u')|  |X;*  sin|a!,„  (u')X.|  -  X;1  sinjaL.  (u')Xk|| 

-  3|X;3sin(ai,s  (u')X.|  -  X;3 cos[a'„„ fu')Ah||>.  for  |u'|  S  W/2.  (Ill 
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For  these  equations,  note  that 


\q  m  Xo/(Ao  —  AX), 

(12) 

\b  *  XqAXo  +  AX), 

(13) 

Si(i)  *  r^dy. 

(14) 

Jo  y 

oJAMu')  *  ir(n  -  n')(l  -  2 u'/W), 

anti(u')  *  tr(n  -  a 'Ml  +  2u’/W), 

0!SJ(u')-»(n  +  t»')(l-2u7W), 

o)in  (u')  *  »(n  +  n')(l  +  2u'/W), 

(15) 

and 

a'„2„>  <u')  -  r[n(l  -  Zu'IW)  -  n'(l  +  2u'/W)|, 
a'*i  (u’)  -  i[n(l  +  2u'/W )  -  n‘ (1  -  2u'/W)|, 
aSU-(u')  -  r[n(l  -  2 u'/W)  +  n'(l  +  2u’/W)\, 
a'*l  (u')  *  »ln(l  +  2u‘/W)  +  n'(l  -  2u'/VF)|.  (16) 

Equations  (10)-(16)  provide  the  mathematical  basis 
for  the  evaluation  of  the  temporal  coherence  require¬ 
ment  of  a  smeared  image  deblurring  process.  Plots  of 
the  normalized  deblurred  image  irradiance,  defined  as 
a  function  of  u'  and  the  bandwidth  AX,  are  shown  in  Fig. 
2.  Xo  **  5461 A  was  used  for  the  calculation.  It  is  evi¬ 
dent  that  if  the  light  source  is  strictly  coherent  (i.e.,  AX 
=»  0),  the  deblurred  point  image  is  infinitesimal.  Also 
note  that  the  degree  of  deblurring  decreases  as  the 
temporal  coherence  of  the  light  source  reduces.  The 
deblurred  length  A IV  represents  the  spread  of  the  de¬ 
blurred  image  irradiance.  It  is  formally  defined  as  the 
separation  between  the  10%  points  of  Z(u').  The  de¬ 
blurred  length  AIV  can  be  shown  to  decrease  mono- 
tonically  with  AX.  Plots  of  AW  as  a  function  of  the 
spectral  bandwidth  AX  for  various  values  of  smeared 
length  W  are  shown  in  Fig.  3.  It  may  be  shown  that  the 
deblurred  length  AW  is  linearly  proportional  to  the 
smear  length  W  for  a  given  value  of  AX.  The  greater 
the  smear,  the  more  difficult  the  deblurring  process.  In 
principle  this  may  be  corrected  by  decreasing  the 
spectral  bandwidth  of  the  light  source.  Table  I  nu¬ 
merically  summarizes  the  preceding  analysis.  The 
value  of  AX  can  be  regarded  as  the  temporal  coherence 
requirement  for  the  deblurring  process. 


B.  Spatial  Coherence  Requirement 
The  relationship  between  the  source  size  and  the 
image  intensity  distribution  is  the  key  factor  in  the 
calculation  of  the  degree  of  spatial  coherence  require¬ 
ment  for  a  particular  image  deblurring  system.  The 
following  analysis  assumes  a  1-D  processor  for  sim¬ 
plicity. 

Assume  the  intensity  distribution  of  the  light  source 
is  given  by 

T(*o)  ■  tect(xo/A*).  (17) 

If  an  extended  monochromatic  light  source  is  used,  the 
output  intensity  distribution  of  Eq.  (1)  becomes 


14 


Fig.  2.  Output  intensity  distribution  of  the  deblurred  image.  AX. 
spectral  bandwidth  of  the  light  source;  W,  smeared  length. 


Spectral  BondwKtm  SXUI 

Fig.3.  Plots  of  the  deblurring  width  A  Was  a  function  of  the  spectral 
bandwidth  of  the  light  source  AX  for  various  values  of  smeared  light 
W. 


V^syA  S  ^ 

A  W/W  1/20  1/15  1/10  1/8  1/5 

AX(A>  270  400  640  750  990 


Hu’)  *  J*  rect(xo/A*)|A(u',xo;X)|Idxo.  (18) 

where 


and  *  denotes  the  convolution  operation. 
Equation  (19)  can  be  reduced  to  the  form  of 
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(  2/A«/  £  (-l)',neapO'2xau7W)  jexpj/  tinc(a  +  n)  -  expj-;  ^-u'jsiiuHa  -  n)|,  for|u'|  >  W/2, 

J)  (-l)"n  *xp(J2rau 7 WO  ( — - — 11  -  exp[/2x(a  -  n)u'/W]  cotx(a  -  n)|} 

I  '  »-i  \a-n  I 


V  + - II  -  exp[/2x(a  +  n)u7W]  cosx(a  +  n)|,  for  |u'|  <  W/ 2, 

a  +  n 

where  a  -  %oW/(\f). 

By  substituting  Eq.  (20)  into  Eq.  (2)  and  performing  and 
the  required  mathematical  evaluation,  the  deblurred  st( 

image  irradiance  is  found  to  be  (a)  for  ju'|  >  W/ 2:  * 


“"V  umivmvumimvw  vTcuwavu/U)  ucuiuiicu  ot/ _  a  \  /n  «/ur  .  t  tl — /  a  v  ,  u 

irradiance  U  found  *  b.  (a)  for  M  >  Stt  ^ 

ig/3\9  (98) 

/•“(u'.At)  m  —22  £  £  |*«>,  (At)  coe(2x(n  -  n')u’IW\  Cf(m,At>  -  C,|t|2u7W  A  1 1  (a(Ai)  +  m|| 

+  ♦«,<*>  ca*(2x(n  +  n')u7W]|.  (21)  '  *  11 1  -S(A,)  +  m|  1 

xt  a  ■  av  ..  ..  -  „  .  Thus  the  spatial  coherence  requirement  of  an  image 

i  eqUatl°n’  f°,l0Wmg  parametere  deblurring  process  may  be  evahited  by  using  Eqs.  (21) 
^  and  (26).  The  output  intensity  distribution  of  this 

Sy(m,n»)-s,  |2x|a(A»)-m)li-sj|2x[-£f(A.)-m||.  (22)  process  is  plotted  in  Fig.  4. 

In  rlHitin  A  r  Recall  the  definition  of  the  deblurred  length  AW 

8  stated  to  be  the  separation  between  the  10%  points  of 

C,(mAi)  -  Cj[2x|a(A»)  +  m|]  -  C,[2x|5(A*)  -  m|],  (23)  the  image  irradiance  I(u').  Figure  5  shows  the  plots  of 


where 


-’■X'  T"' 
CMmS'ir-a' 

S(At)  -  ASW/(2V), 


AW  as  a  function  of  the  source  size  As  and  the  smear 
length  W.  From  this  figure  it  can  be  seen  that  when  the 
spatial  width  of  the  light  source  increases  beyond  a 
critical  size  Ase>  the  deblurred  length  becomes  inde- 


tttwrt-jaLLfclM 

n  -  n  (  in  -  a(At)J 

-  [<-l)"C,(n.At)  -  (-l)"'C/(n',Ai))j  .  (28) 

.  f»n'  (n  +  a(At)| 

♦ffl-(Ai)  “  T7  1  ln  ■  - ]  : 

n  +  n  |  Jn  —  a(At)J 

-[(— l)“Cy(n,Af)  -  <-l)"'Cy(n'.At)|J  • 

(b)  for  |u'|  5  W/2  we  may  similarly  show 


« 

v'.tt) 

......  if 

'''■•if 

Vv/ 

i 

M  U 

i 

\\ 

Z1 

V 

Fig.  4.  Output  intensity  distribution  of  the  deblurred  image  for 
various  values  of  the  source  site  At. 


where 


/(,,(u'.At>  ■  £  £  !♦(£.( At)  +  ♦!S  (A»)  sin[2r(n'  -  n)u7W) 

”  F  t%  n' 

+  *1S  (A»)  cos[2x(n'  -  n)u'/W]  +  +JJ!(At)  cos|2x(n  +  n')u7W]|, 


n3-n'2  In  +  a(A#)J 

+  l-D"  ,.lCt(n'A»)  +  CJ(»',A*)), 
n*  ~  n  * 

+!il(A.)  -  —  |S>(n,At)  +  S»(n,At)], 

n*-  n'2 

♦SS-(As)  -  (- 1 I"*1  -jj-jj— ^ [C» (n At )  +  C»  (n,At)| 

+ ^l|n  [irSSl + [Ci{n’*)  -  c',n'A,)| 
****'  “  rb||n  &TwS)l  “  ,C/(',lA,) + C'(H  • 
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pendent  of  As  and  equal  to  W.  Table  II  provides  a  brief 
numerical  summary  of  the  key  parameters  for  the  de¬ 
termination  of  the  spatial  coherence  in  image  deblurr¬ 
ing.  From  this  table  it  is  evident  that  the  spatial  co¬ 
herence  requirement  is  inversely  proportional  to  the 
smear  length  W.  That  is,  the  longer  the  smearing 
length,  the  smaller  the  source  size  is  required. 

m.  Coherence  Requirement  for  Imago  Subtraction 

Wu  and  Yu18  have  reported  an  image  subtraction 
procedure  employing  an  extended  incoherent  source 
and  a  source  encoded  mask.  Here,  we  shall  discuss  the 
temporal  and  spatial  coherence  requirements  of  the 
image  subtraction  process  and  show  that  they  depend 
on  the  size  and  spectral  bandwidth  of  the  light 
source. 

The  output  intensity  distribution  can  be  derived  by 
calculating  the  propagation  of  the  mutual  coherence 
function  through  an  optical  system  that  utilizes  an  en¬ 
coded  extended  source.  With  reference  to  Fig.  6  the 
encoded  source  intensity  distribution  may  be  described 
by 

iV 

■yUoO'o)  *  r*ctUo/d)  •  £  i(xo  -  nD),  (29) 

where  •  denotes  the  convolution  operation,  D  is  the 
spacing  of  the  2N  +  1  encoding  slits,  and  d  is  the  width 
of  each  slit.  Image  subtraction  is  essentially  a  1-D 
processing  operation  and  will  be  analyzed  as  such. 

The  encoded  mutual  coherence  function  at  the  input 
plane  Pi  can  be  written 

J(u\  -  u»; A) «  J*  _  tUoJ'o)  eipj-/  ^<«i  -  «*>«ojdjto,  (30) 
or,  in  regard  to  Eq.  (29), 

..  d<«  i  ~  uj)  f  .  2v  .  .  J 

J(ui  -  ui;X)  -  £  tine - — - (ui-uj)nD  • 

«— N  V  IV  J 

(31) 

As  shown  in  Fig.  6  two  object  transparencies,  A(u)  and 
B(u),  are  placed  at  the  input  plane  P\.  These  can  be 
represented  by 


lass  —set  at  6pai«l  CWmw  WsejWjwse I 


1/20 

1/15 

1/10 

1/5 

0.5  mm 
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0.38 
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1  mm 

0.1 

0.18 

0.26 

0.40 

2  mm 

0.06 

0.06 

0.12 

0.18 

•  M  M  M  M  l»  I*  14 


Ism*  Silt  MM) 

Pig.  5.  Plot*  of  the  demurring  width  as  a  function  of  the  eouice  size 
for  various  values  of  the  smear  length  W. 


t(u)  -  A(u  -  H)  +  B(u  +  H),  (32) 

where  the  separation  between  transparencies  is  2 H. 
The  mutual  coherence  function  in  the  P2  plane  is  simply 
the  Fourier  transform  of  J0  (ui  ~  U2;X)  convolved  with 
the  spectrum  of  the  amplitude  transmittance  function 
t(u).  The  image  subtraction  process  requires  that  a 
sinusoidal  grating  of  spatial  frequency  uo  be  placed  in 
the  Fourier  plane.  Thus  the  mutual  coherence  function 
behind  this  grating  is  given  by 

«/(*iJtj;X)  ■  J  W(ui  -  Ut;A)«Ui)t*<«2)|(l  +  C  sin2rwoXi) 

X  (1  +  C  sin2vuo*j).  (33) 

A  final  Fourier  transform  operation  due  to  wavelength 
X  will  give  the  intensity  distribution  at  the  output  plane 
P* 


f(u';A)  -  JJ*  ^  JUutiX)  asp  ~  (»i  -  *j)u'  dx\dxt.  (34) 


or 


Fig.  6.  Partially  coherent  optical  proceesing  system  for  image  sub¬ 
traction:  St,  light  source;  Po,  ancodad  extended  source  plans; 
Pi,  input  plana;  PSl  Fourier  plana;  P*  output  ptana;  L,  achromatic 

lenses. 


i 
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/<u';A)  -  N{\A{u'  -  H) I1  +  | B(u'  +  ff)|2]  +  —  |2  sin(2d«o) 

4 

x  R«(A(u'  -  H  +  A/u i0>B(u'  +  H  -  A/«o)l  -  |A(u'  -  H  +  A/<u0)j2 

-  |B(u'  +  H  -  A/u>o)|2  -  |A(u'  -  H  -  A/«0)|2  -  |B(u'  +  H  +  A/«0)|2I.  (36) 


where  Re  denotes  the  real  part  of  |  |.  From  the  above 
equation  we  see  that  there  are  six  diffracted  image  terms 
at  the  output  plane.  If  we  only  consider  the  diffracted 
images  around  the  optical  axis  of  the  output  plane,  from 
Eq.  (35)  we  have 


m  dmA(u')ldu'm, 
B,"°(u')  -  dnB{u')ldu'm. 

we  find  that 


(40) 


f,0,(u';A)  “  |A(u'  -  H  +  A/u>o)|2  -  2  sinc(2du>o) 

x  R«(A(u'  -  H  +  A/a)0)B(u'  +  H  -  \fv0)\ 

+  |B(u'  +  H  -  A/«0)|2.  (36) 

where  u>o  *  H/\f. 

If  the  slit  size  d  equals  zero,  the  analysis  reduces  to 
the  case  of  strictly  spatial  coherence  with  the  intensity 
distribution  given  by 


/,0,(u';A)|j.0  -  (A(u')  -  B(u')]2.  (37) 

Equations  (36)  and  (37)  show  that  a  high-contrast 
subtracted  image  requires  a  strictly  spatial  coherent 
system,  and  that  the  image  quality  will  decrease  as  the 
slit  width  d  increases. 

To  analyze  the  case  of  partial  coherence  assume  the 
light  source  has  a  uniform  spectral  bandwidth  and  that 
the  spectral  response  of  the  detector  is  also  uniform  [i.e., 
s(A)  =  k;  C(A)  =  fe].  Then  the  image  intensity  distri¬ 
bution  at  the  output  plane  »«ay  be  given  by 

/•4V2 

/'°><u')  «  J  ||A(u'  +  AY«o)|2  -  2  «inc(2duo) 

X  RefA(u‘  +  \'fu>o) 

x  B(u'  -  A'/uo)|  +  |B(u'  -  A'to,)|aWA\  (38) 

where  A'  ■  A  -  Ao,  Ao  is  the  center  frequency,  and  AA  is 
the  spectral  bandwidth  of  the  source;  u)0  is  the  spatial 
frequency  of  the  grating.  Equation  (38)  may  be  sim¬ 
plified  by  using  a  Taylor  series  expansion  for  the  input 
object  functions.  Thus  for 


A(u' +  A'/u>o)  *  A(u')  +  £ 

mm  I  ml 

B(u' -  A'/uin)  «  B(u')  +  f)  — -  BiUHu'H\'fwa)m, 

m»i  m! 

where 


(39) 


lot  (b) 

Fit  7.  Apparent  modulation  transfer  function  for  a  partially 
coherent  imace  subtraction:  (a)  basic  frequency;  (b)  second 
harmonic. 


-  [|A(u)|2  -  2  sinc(2du>0)A(u')B(u')  +  |B(u')|2lAA 

-  »inc(2du/o)B,m,(u')] 
x  [A(u')  +  B(u')l(AA)m+* 


mm-  2lm+",-"(m  +  m'+ l)m!m'! 

m  +  m'  •  even 
m  +  0 
m '  *  0 

X  (/m0)*+"  [Al"’,(u')A,",'>(u') 

-  (-I)™'  sinc(2dwo)A ,",,(u')B<l"  )(u') 

+  (-l)",+"''B<"'(u')B'"'»(u')](AA)"+"  ♦>.  (41) 

This  equation  shows  that  a  high-contrast  subtracted 
image  can  be  obtained  with  object  transparencies  of 
moderately  low  spatial  frequency  content.  In  addition 
this  equation  may  be  used  to  compute  the  spectral  re¬ 
quirement  of  the  light  source. 

An  example  analysis  will  be  presented  to  develop  the 
modulation  transfer  function  (MTF)  equations.  These 
will  be  used  to  determine  the  temporal  and  spatial  co¬ 
herence  requirements  of  the  image  subtraction  process. 
Assume  that  the  input  object  transparencies  are  given 
by 

A(u)-1,  B(u)-‘6[1  +  Coco»(2irfc>u)],  (42) 

where  Co  is  the  contrast  of  the  sinusoidal  gTating. 

In  the  case  of  strictly  coherence  (i.e.,  AA  *  d  m  0)  the 
subtracted  image  produces  a  contrast  reversed  image 
with  intensity  distribution 

70(u')  »  lA(u')  -B(u'))2,  (43) 

or 

l  CK  C  CK 

7o(u')  “  -  —  — - co»(2»uiu') - -cob2r(2uj)ii'.  (44) 

4  8  2  8 

For  a  partially  coherent  case  due  to  Eq.  (38),  the  in¬ 
tensity  of  the  subtracted  image  can  be  shown  to  be 

fo(u  )  “  -  ainc(2rfwo)jAA  +  C’ol'fe  -  tinc(2duo)] 

x  »inc(/tax>>oAA)AA  cot(2ruxi') 

Cl 

4  —*inc(2/<iXdoAA)AA  co*2t(2w)u'.  (46) 

O 

Note  the  addition  of  a  second  harmonic  term 
[cos2x(2w)u'|  to  the  basic  frequency.  The  MTF  is 
defined  as  the  ratio  of  the  contrasts  of  the  input  and 
output  sinusoidal  objects.  These  are  given  by 
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_ _  |1  -  2  sinc(2dwo)](2  -  Cj)  sinc(/u<voAA) 

MTF(U” - lO-d-8'.^MW - 


(46) 


MTF(2u>) 


(2  -  Cj)  sin(2/t»xdoAA) 
10  -  Cl  -  8  sinc(2duo) 


(47) 


As  previously  stated  Eqs.  (46)  and  (47)  will  allow  the 
evaluation  of  the  temporal  and  spatial  coherent  re¬ 
quirements  for  image  subtraction. 


A.  Temporal  Coherence  Requirement 

The  case  of  strictly  spatial  coherence  will  be  discussed 
first.  This  requires  that  the  slit  width  d  approach  zero. 
Equations  (46)  and  (47)  will  then  become 

MTF(oi)  ■  sinc(/b>uoAA),  (48) 

MTF(2oi)  -  »inc(2fwwoAA).  (49) 

The  normalized  MTF  curves  of  the  basic  and  harmonic 
frequencies  are  shown  in  Fig.  7.  It  is  obvious  that  the 
contrast  of  the  subtracted  image  decreases  monotoni- 
cally  as  a  function  of  the  object  spatial  frequency. 
However,  the  MTF  of  the  subtracted  image  decreases 
as  the  spectral  bandwidth  of  the  light  source  increases. 
In  other  words  the  quality  of  the  subtracted  image  im¬ 
proves  as  the  spectral  bandwidth  of  light  source  and  the 
spatial  frequency  of  the  object  decrease. 

Let  u)c  be  the  cutoff  spatial  frequency  where  the  MTF 
decreases  to  a  minimum  value  Cm  as  shown  in  Fig.  7. 
The  value  Cm  depends  on  the  maximum  resolution  of 
the  output  recording  material.  Figure  8  shows  the 
functional  relationship  of  the  cutoff  frequency  «t.  and 
the  spectral  width  AA  for  various  values  of  Cm.  It  is 
possible  to  determine  the  spectral  bandwidth  require¬ 
ment  AX  from  this  figure.  The  relationship  between 
the  cutoff  frequency  wc,  the  spectral  width  AA,  and  the 
separation  between  two  input  transparencies  H  is  shown 
in  Fig.  9.  Note  that  the  spectral  bandwidth  required 
for  a  given  cutoff  frequency  decreases  with  increasing 
separation  H.  Table  III  illustrates  the  dependence  of 
AA  on  <*>c  and  H.  The  focal  length  of  the  Fourier 
transform  lens  selected  is  /  ■  300  mm  for  calculation. 
It  is  clear  from  the  table  that,  as  the  spatial  frequency 
and  object  separation  increase,  the  spectral  bandwidth 
of  the  light  source  must  decrease. 


Fig.  8.  Relationship  between  the  cutoff  frequency  wo  and  the  spec¬ 
tral  bandwidth  of  the  light  source  AA  for  different  minimum  desirable 
contrasts  C„. 


Fig.  9.  Relationship  between  the  cutoff  frequency  end  the  spectral 
bandwidth  of  the  light  source  AA  for  various  values  of  separation  H . 
2 H  is  the  main  separation  between  the  input  object  transparencies. 
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B.  Spatial  Coherence  Requirement 

Consider  the  case  of  perfect  temporal  coherence  (i.e., 
AA  *  0)  and  partially  spatial  coherence,  where  Eqs.  (46) 
and  (47)  are  of  the  form 


MTF, 


MTF2 


(2  -  Cj)(l  -  2  sinc(2d<oo)) 
10  -  Cl  -  8  sinc(2<fwo) 

2-q 

10  —  Cl  -  8  sinc(2d<vo) 


(50) 

(51) 


Note  that  the  MTFs  are  independent  of  the  object’s 
spatial  bandwidth.  The  above  equations  are,  however, 
dependent  on  the  slit  width  d.  This  requires  that  the 
grating  be  precisely  designed  to  match  the  separation 
of  the  object  transparencies,  i.e., 


<*>  -  H/IX/o).  (52) 

The  plots  of  the  MTF  vs  the  separation  H  for  Eqs.  (50) 


r%.  10.  Apparent  modulation  transfer  function  vs  the  separation 
H  for  different  desirable  contrasts  0„ 
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and  (51)  are  shown  in  Fig.  10.  It  is  obvious  that  to  ob¬ 
tain  a  high-contrast  subtraction  image  the  separation 
H  must  be  reduced.  However,  decreasing  the  separa¬ 
tion  between  the  object  transparencies  limits  the  size 
of  input  objects  to  be  processed.  The  relationship  be¬ 
tween  the  MTF,  the  object  transparency  separation  H, 
and  the  slit  width  d  is  numerically  presented  in  Table 
IV,  where  the  focal  length  of  the  Fourier  transform  lens 
was  assumed  to  be  300  mm. 

Note  that  this  table  indicates  the  necessity  of  a  very 
narrow  source  size  to  achieve  an  adequate  MTF. 
However,  to  obtain  a  high-intensity  narrow  source  size 
is  difficult  to  achieve  in  practice.  The  problem  can  be 
solved  considerably  if  source  encoding  techniques  for 
image  subtracting  are  used.17-18 

A  multislit  source  encoding  mask  is  used  for  this  il¬ 
lustration.  The  spatial  period  of  the  encoding  mask 
should  be  precisely  equal  to  that  of  the  diffraction 
grating  G  (i.e.,  D  ■  l/o>0).  The  spatial  coherence  re¬ 
quirement,  although  independent  of  the  slit  size,  is 
governed  by  the  ratio  of  the  slit  width  to  the  spatial 
period  of  the  encoding  mask,  i.e.,  d/D.  The  ratio  d/D 
must  be  relatively  small  to  achieve  a  high  degree  of 
spatial  coherence.  The  dependence  of  the  MTF  on  d/D 
is  shown  in  Fig.  11.  It  is  obvious  that  the  subtraction 
effect  ceases  when  the  MTF  approaches  zero,  i.e.,  d/D 
«  0.3.  A  few  numerical  examples  are  presented  in 
Table  V. 

IV.  Summary 

We  have  shown  that  the  temporal  and  spatial  co¬ 
herence  requirements  for  some  partially  coherent  op¬ 
tical  processing  operations,  namely,  image  deblurring 
and  image  subtraction,  can  be  determined  in  terms  of 
the  output  intensity  distribution.  For  image  deblurring 
the  temporal  coherence  requirement  depends  on  the 
ratio  of  the  deblurring  width  to  the  smeared  length  of 
the  blurred  image.  To  obtain  a  higher  degree  of  de¬ 
blurring  a  narrower  spectral  width  of  the  light  source 
is  required.  For  example,  if  the  deblurring  ratio  Aw/w 
is  0.1,  the  spectral  width,  AX,  should  be  <640  A. 

For  the  spatial  coherence  requirement  the  image 
deblurring  depends  on  both  the  deblurring  ratio  Aw/w 
and  the  smeared  length  w.  If  the  deblurring  ratio  Aw/w 
*  1/10  and  w  *  1  mm,  a  slit  source  <0.26  mm  should  be 
used.  For  a  smeared  image  deblurring  operation  the 
constraints  of  the  temporal  and  spatial  coherence  re¬ 
quirements  are  not  critical,  which  can  be  achieved  in 
practice. 

For  image  subtraction,  the  temporal  coherence  re¬ 
quirement  is  determined  by  the  highest  spatial  fre¬ 
quency  and  the  separation  of  the  input  object  trans¬ 
parencies.  If  the  separation  and  spatial  frequency  of 
the  input  transparencies  are  high,  a  narrower  spatial 
bandwidth  of  the  light  source  is  required. 

For  spatial  coherence  requirement,  the  modulation 
transfer  function,  which  determines  the  contrast  of  the 
subtracted  image,  depends  on  the  ratio  of  the  slit  width 
to  the  spatial  period  of  the  encoding  mask,  i.e.,  d/D.  If 
the  ratio  d/D  is  low,  a  higher  contrast  subtracted  image 
can  be  obtained.  For  example,  with  d/D  ■  0.05,  a  rel- 


TaMt  IV.  Sewed  Six*  for  knap*  SsWracSMi  Under  OMaraal  MTS  and 
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Pig.  11.  Relationship  between  MTF  (at),  MTF  (2w),  and  the  ratio 
of  the  alit  width  to  apatial  period  d/D. 
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atively  higher  MTF  “  0.85  can  be  obtained.  Compared 
with  the  image  deblurring  operations,  the  coherence 
requirements  are  more  stringent  for  the  subtraction 
process. 

Finally,  the  solution  to  the  coherence  requirement  for 
partially  coherent  processing  is  not  restricted  to  the 
application  of  the  deblurring  and  subtraction  operation, 
but  may  also  be  applied  to  any  other  optical  processing 
operation. 

We  wish  to  acknowledge  the  support  by  the  U.S.  Air 
Force  Office  of  Scientific  Research  grant  AFOSR-81- 
0148. 


Appendix 

To  prove  the  following  relationship: 

PS5 1  (-1)n 


Wx 

:  w 


X  tin 


IV  „T, 
2«tiXo 


XIV 


u'|  agn(u 


<A1> 


we  let  x/\f  m  fx  and  X/Xo W  *  l. 

The  problem  then  reduces  to  finding  the  solution 
for 
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v/  f  ■  7~rr  »»P0  2t/xii  'Mix.  (A2) 

J  sui(t(/x) 

Since  there  are  an  infinite  number  of  poles  in  the  real 
domain,  i.e., 

fx  *  nil,  n  *  1,2 .  (A3) 

Equation  (Al)  will  be  evaluated  for  the  cases  of  u'  >  0 
and  u'  £  0.  The  contour  integration,  which  is  taken 
over  the  upper  half  of  the  complex  plane  (see  Fig.  12), 
is  given  by 


IMUI 


|X'.  -  LC  *  XI  «*-'■ !"“w-  -  «■ 

(A4) 

where  x  is  the  radius  of  the  small  half -circles  around  the 
poles,  and  R  is  the  radius  of  the  larger  contour  half¬ 
circle.  Denote 

z  •  R  exp (jf»,  dz  -  jBR  exp(j#)d9.  (A5) 

The  last  term  of  Eq.  (A4)  may  then  be  written 

f  -•  TV,-;  «*p02 *fxu'\dfx  -  f -~T,  exp(j2rzu')dz 
Jr  sin(a7/x)  J  sin(irfc) 

J*  exp (jirRtf  cos#)  exp(-2rR  si n8u‘)B  exp(2jB)  ^ 

R  exp[yir(/t(co8#  +  j  sin#))  -  exp[-;>/fl(cos#  +  j  sin#)] 


lim  f  7  -  np(j2rfxu')dfx  «  0.  (A7) 

*— •  Jr  sin irlfx) 

But  from  Eq.  (A4)  it  is  easily  shown  that  as  R  ap¬ 
proaches  infinity 

f  T~~  ■  txp(j2rfxu')dfx  *■  -2r;  £  K„,  (A8) 

•/—  sin(x7/x)  ««-« 

where 

Rn  .  limfx-  S|  ”  «pO;2?^VO  .  (A9) 

i— f  V  U  sin(rli)  rl1 

Therefore, 

f  :  exp02tr/xu’)d/x  ■  ^  £  (-lVn  sin — u', 

J—  im(tlfx)  i*n.  i  l 

to ru’>0.  (A10) 

The  contour  integral  taken  over  the  lower  half  of  the 
plane  as  shown  in  Fig.  13  is  similarly  given  by 

f  --77—  exp(/2»/x«'W/x  -  -7  £  (~l)"n  sin  — p u'. 
J— sin(*f/x>  (*«- 1  ( 

for  u'  <  0.  (All) 

The  results  of  Eqs.  (A10)  and  (All)  thus  prove  that  the 
equality  of  Eq.  (Al)  holds  true. 
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Abstract.  In  this  paper,  the  general  formulations  of  the  apparent  transfer  function  for  the 
partially  coherent  optical  processor  will  be  derived.  Although  these  formulas  show  that  the 
apparent  transfer  function  is  dependent  upon  the  degree  of  spatial  and  temporal  coherence, 
there  is  actually  more  variability  in  the  spatial  coherence.  We  note  that  the  obtained 
formulas  may  also  be  used  as  a  criterion  in  the  selection  of  source  size  and  spectral 
bandwidth  of  an  incoherent  light  source.  Thus  a  specific  optical  information  processing 
operation  can  be  carried  out  with  an  incoherent  source. 

PACS:  42.30,  42.80 


The  description  of  a  transfer-function  for  a  linear 
spatially  invariant  optical  system  is  an  important 
concept  in  image  evaluation.  The  techniques  of  using 
amplitude  and  intensely  sinusoidal  objects  as  input 
signals  to  determine  the  transfer  characteristic  of  a 
coherent  and  an  incoherent  optical  system  has  been 
investigated  previously  [1, 2].  In  the  past,  the  concept 
of  system  transfer  function  has  been  used  as  a  criterion 
to  evaluate  the  image  quality  of  an  optical  system. 
However,  a  strict  coherent  or  incoherent  optical  field  is 
a  mathematical  idealization.  An  optical  field  that 
occurs  in  practice  consists  of  a  very  limited  degree  of 
coherence,  because  the  electromagnetic  radiation  from 
a  real  physical  source  is  never  strictly  monochromatic 
In  reality,  a  physical  source  cannot  be  a  point,  but 
rather  a  finite  extension  which  consists  of  many  ele¬ 
mentary  radiators. 

The  optical  system  under  the  partially  coherent  regime 
has  been  studied  by  Becherer  and  Parrent  [3],  and 
Swing  and  Clay  [4].  They  have  shown  that  there  are 
difficulties  in  applying  the  linear  system  theory  to  the 
evaluation  of  imagery  at  high  spatial  frequencies. 
These  difficulties  are  primarily  due  to  the  inapplica¬ 
bility  of  the  linear  system  theory  under  partially 
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coherence  regime.  Nevertheless,  using  a  sinusoidal 
analysis,  an  apparent  transfer  function  for  a  partially 
coherent  optical  system  was  obtained.  The  result  is 
appreciably  different  from  those  that  would  have  been 
obtained  from  a  linear  system  concept,  either  in  in¬ 
tensity  or  in  complex  amplitude.  In  a  more  recent 
paper,  Dutta  and  Goodman  [5]  described  a  procedure 
for  sampling  the  mutual  intensity  function  so  that  the 
image  of  a  partially  coherent  object  can  be 
reconstructed. 

In  this  paper,  we  shall  show  that  an  apparent  transfer 
function  for  an  optical  information  processing  system 
can  be  derived  from  a  partially  coherent  illuminator. 
We  shall  show  that  the  apparent  transfer  function  of  a 
partially  coherent  optical  processor  is  dependent  upon 
the  temporal  and  spatial  coherence  of  the  light  source. 
We  shall  also  show  that,  the  concept  of  transfer 
function  is  a  valuable  one  that  can  be  used  as  a 
criterion  for  selecting  an  appropriate  incoherent  light 
source  for  a  specific  information  processing  operation. 


1.  Apparent  Transfer  Function  foe  Temporal  Coherence 

In  a  recent  paper  [6],  we  have  evaluated  the  coherence 
requirement  for  a  partially  coherent  optical  infor¬ 
mation  processing  system.  We  have  shown  that  the 
temporal  and  spatial  coherence  requirements  are,  re- 


0721-7269/82/0028/0359/S01.60 


360 


23 


S.  L.  Zhuang  and  F,  T.  S  Yu 


Fig.  1.  A  partially  coherent  optical  information 
processor.  (L:  achromatic  transform  Ians.  p0. 
source  plane:  p, :  input  plane:  p2:  spatial  fre¬ 
quency  plane:  ps :  output  plane) 


spectively.  dependent  upon  the  spectral  bandwidth  and 
the  source  size  of  the  light  source.  In  the  analysis  of  the 
temporal  coherence  requirement,  we  have  assumed  an 
infinitely  small  light  source  with  finite  spectral  width. 
And  for  the  evaluation  of  the  spatial  coherence  re¬ 
quirement.  we  let  monochromatic  source  with  finite 
extent.  In  the  evaluation  of  the  apparent  transfer  for  a 
partially  coherent  processing  system,  we  shall  also  use 
these  basic  approaches. 

We  shall  now  evaluate  the  relationship  between  the 
apparent  transfer  function  of  an  optical  system  and  the 
spectral  bandwidth  of  a  light  source.  Let  us  now  refer 
to  the  partially  coherent  optical  information  system  of 
Fig.  1.  We  note  that  a  diffraction  grating  is  used  at  the 
input  signal  plane  (u,v).  The  purpose  of  using  a 
diffraction  grating  at  the  input  plane  is  to  disperse  the 
input  signal  in  the  spatial  frequency  plane  (x,  y)  so  that 
a  high  spatial  coherence  can  be  achieved  in  the  Fourier 
domain  [7],  Thus  the  input  signal  can  be  processed  in 
a  complex  amplitude,  rather  than  in  intensity,  with  a 
broad-band  incoherent  or  white-light  source.  By  using 
the  partial  coherent  theory  [8],  we  have  computed  a 
general  output  intensity  distribution  for  the  partially 
coherent  processor,  as  shown  in  the  following  [7] : 


Ao+JA/2  x 

Hu',v')=  j  fj  y(x0,.v0)s<A)cU) 

Ao  “  JA,  2  -  x 

jj  T(x0  +  x  -  A/p0,>’o+y)/(x,>’)e  V"  >v)dxdy\ 


■ dx0dy0dX , 


(1) 


ing,  can  be  written  as: 

f(u)  =  (l  +ccos2nwuMl  +cos27tp0u).  (2) 

where  a>  is  the  spatial  frequency  of  sinusoidal  signal 
and  c  is  the  contrast.  For  simplicity,  we  use  a  one 
dimension  analysis.  Thus,  the  output  intensity  distri¬ 
bution  of  Eq.  (1)  can  be  reduced  to: 

Ao  +  4A/2  x 

I(u',AX)  =  J  {  y(x0)|/t(x0,u';  A)|2dx0dx,  (3) 

Ao  -  JA/2  —  x 

where 

x  _ .2*  . 

/4(x0,  u  ;  /.)  =  J  T(x0,x:  A)/(x)e  >if*“dx.  (4) 

-  X 

is  a  complex  amplitude  function. 

In  order  to  achieve  a  higher  temporal  coherence 
requirement  for  the  complex  filtering  process,  we  limit 
the  spatial  filter  /(x,y)  to  a  finite  extension  in  the  x 
direction,  i.e., 

/(x)  =  rect|*  ^-°j,  (5) 

where  Ax  is  the  spatial  width  of  the  filter,  and  the  filter 
is  centered  at  x  =  Xfp0  in  the  spatial  frequency  plane. 
Since  the  Fourier  spectrum  of  the  input  signal  at  tH? 
spatial  frequency  plane  p2  is: 

X 

T(x0,x;x)=  |  (1 +ccos2ji«u) 

-®  ,z«  +  , 

(1 +cos2np0u)e  >if  *°  *'du  (6) 


where  A0  and  AX  are  the  central  wavelength  and 
spectral  width  of  the  light  source,  y(x0,  y0)  is  the  source 
intensity  distribution,  s(X)  and  c(x)  are  the  source 
spectral  distribution  and  the  spectral  sensitivity  of  the 
recording  materia],  respectively.  /(x,y)  is  the  filter 
function  at  the  Fourier  plane,  T(x,y)  is  the  Fourier 
spectrum  of  the  input  signal  and  p0  is  the  spatial 
frequency  of  this  diffraction  grating.  For  simplicity,  we 
assume  that  both  s(A)  and  cU)  are  constants  which  can 
be  ignored  in  (1).  We  now  use  a  sinusoidal  input  object 
transparency  to  evaluate  the  apparent  transfer  func¬ 
tion  of  the  optical  processing  system.  The  overall  input 
signal  transmittance,  that  includes  the  diffraction  grat- 


we  obtain : 


lxn  +  x  \ 

1  . 

fx0+x  \ 

O 

Ci. 

1 

> 

+  2* 

V  +Po) 

♦i'ftr--**) 
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(7) 


where  p0  5>  <o. 

In  optical  signal  processing  we  focus  only  the  first- 
order  terms  of  (7)  in  the  following 

r'ftf  + 


+j4' 


(x0  +  X  _ 

l  4 


(Po+w) 


(8) 


By  substituting  (5)  and  (8)  into  (4),  we  have: 
A(x0,  u‘;k)  =  ]  rect  |—  >k4u  dx 

•IW^+ 

*JL— 


2« 

^“dx 


where  •  denotes  the  convolution  operation. 

We  now  consider  the  effect  due  to  temporal  coherence. 
Let  us  now  assume  that  the  incoherent  source  is  a 
point  source,  i.e.,  y(x0)*<5(x0),  thus  (3)  becomes 

10+42/2 

/(u',dA)=  J  /t(u';  *)/«•(«' U)d2 

io-ai/2 

2o+  42/2  1 1  „2 

J 


{£} 
my- 


2o-42/2 


+  -cos(2jttuu’)rect 


+  —  cos(2n(2w)u')  rect 

O 


(10) 


where  Ak  is  the  spatial  bandwidth  of  the  light  source. 
In  the  following,  we  shall  attempt  to  evaluate  the 
output  intensity  distribution  of  (10)  in  the  following 
separated  cases : 

/i)  for  o><  — r  —  - — ,  the  output  intensity  distri- 
(2/.0  +  J/.)f 

bution  of  Eq.  (10)  becomes 
I"\u.Ak) 


1  c2  c  c 2 

-  +  —  +  -cos(2!Kok')  +  cos(2*2«u') 

4  8  2  8 


dx,  (11) 


where  the  image  contrast  for  both  basic  and  second 
harmonic  frequencies  are 


vVW)= 


and 


>'"(«>)= 


(iij  for 


2  +  c2  ’ 


2  +  e2  ’ 
Ax 


(12a) 


(12b) 


<(!)< 


Ax 


f{2k0  +  Ak)  fako-Ak) 
2 


— ,  we  obtain 


^u,Ak)=lAk+^-k0+f) 

+  C-  cos(2«tw')(^ 

+  cjCOs(2n2wu')(^)-^fy  (13) 


and 


c  j2«(^+<4« 

+  -e  W  ^  rect 

fx0  +  kfwA 

y\2lM = 

4 

l  Ax  \ 

C  i2"(77 -«)» 

+  -e  ’  rect' 

fx0-2>l\  2,«.’ 

4 

1  dx  J; 

//'(«)= 

(9) 

(Hi)  for  o>> 


4 c{Ax  —  2fuikg  +  fa>Ak) 
4fwAk  +  c2(dx  —  2/eoAo +fa)Ak)  ’ 

c2(  Jx  -  2fwk0  4-  foiAk) 
4/tuAk  +  c2{Ax  -  2fa>k0  +  fcoAk)  * 

Ax 


A2ko  -Ak) 


,  (10)  becomes 


rV.dA)-  \Ak. 

4 

and 


(14a) 


(14b) 


(15) 


(16) 


We  also  note  that  for  the  case  of  monochromatic  point 
source,  the  output  irradiance  distribution  is 


c2  c2 

/(m)*  1  +  —  +2ccos(2jkou)+  —  cos(2tt2tou) , 


(17) 
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Fig.  2.  Temporal  coherence  apparent  transfer  function  as  a  function 
of  input  signal  frequency  a>,  for  various  values  of  spectral  width  AX 


MTF(W) 


Fig.  3.  Temporal  coherence  apparent  transfer  function  as  a  function 
of  to,  for  various  of  spatial  width  of  filter  Ax 


where 

yo»-  40 

7l  2  +  cJ  ’ 

(18a) 

and 

yo>_  c2 

72  2  +  c1 

(18b) 

Because  the  apparent  transfer  function  can  be  obtained 
by  the  ratio  of  the  contrast  of  the  output  basic 
frequency  signal  to  the  contrast  of  the  input  signal,  we 
have 

11  for  a< 


Which  apparent  transfer  function  reduces  to  a  function 
of  p0  and  Ak. 

We  now  illustrate  the  dependency  of  the  apparent 
transfer  function  upon  the  spectral  width  of  the  light 
source  Ak,  the  spatial  width  of  the  filter  Ax,  and  the 
spatial  frequency  of  the  grating  p0.  Figure  2  shows  the 
plots  of  the  apparent  transfer  function  (MTF)  as  a 
function  of  the  input  signal  frequency  w  for  various 
values  of  spectral  bandwidth  Ak.  From  this  figure  we 


Ax 

t0  +  Ak) 


MTF(a>)  = 


(2  +  c2)(dx  —  2f cak0  +  f coAk) 
4fa)Ak+c1(Ax  —  2fa)k0  +  fa)Ak) 


for 


Ax 


f(2k0  +  Ak) 


Ax 


f(2k0-Ak) 


IP  for  w  2; 

which  is  depending  upon  the  spatial  frequency  of  the 
diffraction  grating  p0,  the  spatial  width  of  the  filter  Ax, 
and  the  spectral  bandwidth  of  the  light  source  Ak. 

We  note  that  the  width  of  the  Alter  Ax  should  be 
chosen  to  coincide  with  the  product  of  the  spatial 
frequency  p0  of  the  diffraction  grating,  the  focal  length 
of  the  achromatic  transform  lens,  and  the  spectral 
width  Ak  of  the  light  source,  Le., 

Ax  =  PofAk.  (20) 


Ax 

m  k0-Ak) 


(19) 


notice  that  the  MTF  is  not  appreciably  affected  by  the 
spectral  bandwidth  Ak  of  the  light  source,  except  for 
some  slight  changes  in  frequency  response.  For  exam¬ 
ple,  an  Ak  becomes  broader,  a  slight  reduction  in  lower 
frequency  response  is  expected.  However,  the  system 
bandwidth  is  somewhat  slightly  broader.  Figure  3 
shows  the  variation  of  MTF  as  a  function  of  w  for 
various  values  of  spatial  width  Ax  of  the  filter  and  a 


By  substituting  this  relationship  of  (20)  into  (19),  we 
have 


MFT(co)» 


for  w< 


(2  -I-  c1)(p0Ak-2o)k0+  (oAk) 


4 <oAk  +  c2(p0Ak—  2wk0  +  a)Ak) 


for 


2k0  +  Ak 
Po * 


for  cu£ 


2A„-I-  Ak 
p0k 


£a>< 


2ka  -  Ak ' 


Po* 

2^ -Ak 


10 


(21) 
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MTF(Ctf) 


Fig.  4.  Temporal  coherence  apparent  transfer  function  as  a  function 
of  co,  for  various  values  of  grating  spatial  frequency  p0 


given  dA.  From  this  figure  we  see  that  the  system 
bandwidth  is  linearly  related  to  that  of  the  spatial 
width  of  the  filter.  In  other  words,  the  larger  the  spatial 
width  of  the  filter  used,  the  wider  the  system  band¬ 
width  which  may  result.  Figure  4  also  shows  the 
dependent  of  the  MTF  upon  the  spatial  frequency  of 
the  grating  p0.  Again,  we  see  that  the  system  band¬ 
width  is  linearly  proportional  to  p0,  as  expected  from 
the  relationship  of  (20).  In  other  words,  a  higher 
frequency  grating  has  the  advantage  of  achieving  finer 
image  resolution.  However,  this  advantage  of  using  a 
higher  spatial  frequency  grating  is  somewhat  com¬ 
pensated  with  the  use  of  a  larger  achromatic  transform 
lens,  to  which  is  generally  more  expensive. 


2.  Apparent  Transfer  Function  for  Spatial  Coherence 

Now  we  shall  determine  the  dependent  of  the  apparent 
transfer  function  upon  the  source  size  As.  To  simplify 
our  analysis,  we  assume  that  the  light  source  is  mono¬ 
chromatic  but  with  finite  extent.  The  output  intensity 
distribution  of  (3)  can  therefore  be  written  as : 

/(u',ds)=  J  y(x0)|d(u',x0)|2dx0.  (22) 


For  simplicity,  we  assume  that  a  uniform  irradiance 
distribution  of  the  extended  light  source  exists,  Le., 


y(x0)«rect  j-^j, 


26 
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where  As  is  the  extended  source  size  with  reference  to 
(4),  we  see  that 

A{u,  Xq)A  *{u',  Xq)  —  ^rect  j^j 

+  2C0s(2ttwu')(rect  1*° 

4  \  ldx-/i0/wj 

[*o-iA0/a)n 

Mx-A0/o>j/ 


+  rect 


+  —  cos(2n2a>u')  rect  \ 

O 


(24) 


J _ _ l 

Idx— 2A0/tt>J 
where  we  have  used  the  following  relationship 

(25) 

By  properly  shifting  the  coordinate  axis,  (24)  can  be 
written  as 


1  i  X  1 

A(u\  x0)A*(u',x0)=  -rect 

+bM 


rr-^-)+rect| 


dx  +  2A0/<a) 


J) 


+  -cos(2jrwu')^rect 


|M 

IdxJ 


+  rect 


(dx  — 2Ap/a>J 
*  -}) 


(- 


+  —  cos(2n2a>u')  rect 

O 


[Ax-lX^fw 
|dx  — 2Ap/(i>}’ 


(26) 


By  substituting  (23)  and  (26)  into  (22),  we  obtained  the 
following  output  intensity  distribution  and  the  basic 
and  second  harmonic  frequencies 
(i)  for  0<ds<dx  —  2Ap/w,  we  have: 


/U)(u',  ds)  = 
A 


+  —  cos(2jr2aw') 

O 


+  y)  +  ^cos(2ita>u') 
As, 


and 

/i“(a>)= 

y(i2*M= 


4c 

2  +  7 


2+c 2 ' 

(ii)  for  Ax  —  2Ap/ci><ds<dx,  we  obtain 
fV.  ds)= 


(27) 

(28a) 

(28b) 


(\+r^s+r6iAx-2^f(a) 


+  -{As  +  dx-2Ap/a>)cos(2jKuu') 


+  —{Ax  —  2Ap/ w)cos(27t2cau') , 


(23) 


(29) 
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J6S 


and 


'/j4'(w)~ 


c1 2(J.x-2/.0/<u) 
(2  +  c2)d.x 


(34b) 


( tv  for  d.x  <  a>A0f  we  again  obtain 

/,5,(h.  Js)  =  constant .  (35) 

and 

•/,5W)"-//W)-0.  (36) 

Notice  that  where  this  cut-off  frequency  of  the  optical 

system  is  to =  Ax!An  f.  Thus  the  apparent  transfer 

function  of  the  optical  processing  system,  for  both  the 
basic  and  second  harmonic  frequencies  are 


the  source  size  increases.  Unlike  the  first  harmonic 
MTF  of  Fig.  5a,  the  cut-off  frequency  tends  to  stay  at 
the  same  values,  although  the  higher  frequency  re¬ 
sponse  decreases  rather  rapidly. 

Figures  6a  and  6b  show  the  basic  and  the  second- 
harmonic  MTF  as  a  function  of  input  spatial  frequency 
a>,  for  various  sizes  of  spatial  filters  Ax,  for  a  given  As. 
From  these  two  figures,  again  we  see  that  the  transfer 
system  bandwidth  is  linearly  related  to  the  size  of  the 
spatial  filter.  Thus  MTF  is  obviously  limited  by  the 
filler  bandwidth  of  the  optical  processor.  The  increase 
of  filter  bandwidth  also  causes  a  reduction  in  temporal 
coherence,  which  in  turn  reduces  the  processing  capa¬ 
bility.  Nevertheless,  an  optimum  processing  capability 


MTF(oi) 


1 

(2  +  c2)(  As  +  Ax-  2  ip/ta) 
4ds  +  c2(4s  +  Ax—  2 a0/<u) 
2(2  +  c2)(Ax-A0fw) 

4  Ax  +  c2{As  4 Ax  —  2A0f ai) 
Ax  —  X0fu) 

Ax 

.0 


for  0  £  As  <  Ax  -  2A0f  a> 
for  Ax  —  2A0fa>&As<Ax 

for  Ax^As%Ax  +  2A0fo) 

for  Ax  +  2X0/oj^  As 
for  Ax£atA0f 


and  similarly : 


MTF(2<o)» 


1 

2(2  +  c2)(Ax-2A0fco) 

4As+c2(As  +  Ax  —  2A0f to) 
2(2  +  c2)(dx-2A0/<a) 

4  Ax  +  c2(ds  +  Ax  —  2  A0f to) 
Ax-2A0foi 
lx 

0 


for  0£As<Ax  —  2A0fw 
for  Ax-  2A0j<u  £As<Ax 

for  Ax£  As<Ax  +  2A0fu> 

for  Ax  +  2A0fa)£  As 
for  Ax^Aq/w. 


(37) 


(38) 


From  the  above  equations,  one  may  see  that  MTF, 
either  for  the  fundamental  or  for  the  second-harmonic 
frequency,  decreases  rather  rapidly  as  the  source  size 
As  increases.  Figure  5a  shows  MTF  as  a  function  of 
input  spatial  frequency  to,  for  various  values  of  source 
size  ds.  From  this  figure,  we  discover  that  the  fre¬ 
quency  response  (i.e.,  MTF)  decreases  quite  rapidly  as 
the  source  size  increases.  In  other  words,  for  a  fixed 
filter  size,  the  smaller  the  source  size  used,  the  better 
the  system  frequency  response.  We  also  see  that,  where 
the  source  size  is  adequately  small  (e.g.,  As  $0.02  mm) 
the  MTF  approaches  the  strict  coherent  MTF,  on  the 
other  hand,  if  the  source  size  is  significantly  large  (e.g., 
As  £  4  mm)  the  MTF  approaches  that  of  the  incoherent 
case.  Figure  5b  shows  the  second-harmonic  MTF  as  a 
function  of  input  spatial  frequency  a>.  From  the  figure, 
we  see  that  frequency  response  decreases  even  faster  as 


can  be  obtained  with  the  appropriate  MTF  for  certain 
optical  processing  operations. 


The  nonlinear  behavior  of  the  partially  coherent  opti¬ 
cal  processor,  when  considering  either  intensity  or 
amplitude  distribution  input  signals,  necessitates  the 
use  of  the  apparent  transfer  function  to  accurately 
predict  the  system  response.  We  have  derived  the 
general  formulas  for  MTF  in  terms  of  the  theory  of 
partially  coherent  light.  These  derivations  indicate  the 
dependence  of  MTF  upon  the  degree  of  spatial  coher¬ 
ence  (Le.,  the  source  size)  as  well  as  the  degree  of 
temporal  coherence  (Le.,  the  source  spectral  band¬ 
width).  MTF  has  been  shown  to  be  less  dependent 
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upon  the  spatial  coherence  requirement  as  compared 
to  its  relationship  with  the  temporal  coherence 
requirement. 

It  has  been  noted  that  the  spatial  bandwidth  of  our 
optical  processor  is  primarily  dependent  upon  the  size 
of  the  filter  Ax,  where  the  filter  is  placed  in  the  Fourier 
plane.  The  transfer  systems  bandwidth  may  be  in¬ 
creased  by  using  a  larger  spatial  filter  Ax.  However,  the 
size  of  the  filter  is  selected  such  that  Ax = p0f  AX,  which 
is  linearly  related  to  the  spatial  frequency  of  grating, 
the  focal  length  of  the  transform  lens  and  the  spectral 
width  of  the  light  source.  A  narrow  spectral  band  AX  is 
necessary  for  most  partially  coherent  optical  infor¬ 
mation  processing  operations.  In  order  to  achieve  the 
required  AX  for  a  wide  strip  of  spatial  filter  Ax  in  the 
spatial  frequency  plane,  a  diffraction  grating  of  suf¬ 
ficiently  high  frequency  p0  at  the  input  plane  is  needed. 
For  example,  for  partially  coherent  processing  with  a 
white  light  source,  a  set  of  narrow  spectral  band  filters, 
each  with  a  spectral  bandwidth  AX,  can  be  used  in  the 
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spatial  frequency  plane.  Finally,  we  stress  that  the 
apparent  transfer  function  which  we  have  obtained  is 
rather  general  and  may  be  applied  to  any  partially 
coherent  optical  processing  system. 
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Abstract 

A  relation  between  spatial  coherence  and 
source  encoding  intensity  distribution  is  pre¬ 
sented.  Since  the  spatial  coherence  require¬ 
ment  is  determined  by  the  signel  processing 
operation,  a  strict  coherence  source  may  not  be 
needed  for  the  processing.  The  advantage  of  the 
source  encoding  is  to  relax  the  constraints  of 
strict  spatial  coherence  requirement  so  that 
the  signal  processing  operation  can  be  carried 
out  with  an  extended  incoherent  source.  The 
effects  of  signal  sampling  on  coherence  require¬ 
ment  is  discussed.  The  advantage  of  signal 
sampling  is  to  Improve  the  degree  of  temporal 
coherence  In  the  Fourier  plane  so  that  the  sig¬ 
nal  can  be  processed  in  complex  amplitude  with 
a  broadband  white-light  source.  Finally, 
experimental  demonstrations  are  Included. 


Introduction 

Since  the  invention  of  laser  (i.e.,  a 
strong  coherent  source^,  laser  has  become  a 
fashionable  tool  for  many  scientific  applica¬ 
tions  particularly  as  applied  to  coherent  optical 
signal  processing.  However  coherent  optical 
signal  processing  systems  are  plagued  with  co¬ 
herent  noises,  which  frequently  limit  their 
processing  capability.  As  noted  by  the  late 
Gabor,  the  Nobel  prise  winner  in  physics  in  1970 
for  his  invention  of  holography,  the  coherent 
noise  is  the  number  one  enemy  of  the  Modern 
Optical  Signal  Processing1.  Aside  the  coherent 
noise,  the  coherent  sources  are  usually  expan¬ 
sive,  and  the  coherent  processing  environments 
are  very  stringent.  For  example,  heavy  optical 
benches  and  dust  free  environments  are  generally 
required. 

Recently,  we  have  looked  at  the  optical  pro¬ 
cessing  from  a  different  standpoint,  A  question 
arleea,  is  it  necessarily  true  that  all  optical 
signel  processing  required  a  coherent  sourest 
The  answer  to  this  question  is  that  there  are 
many  optical  signal  processings  that  can  be 
carried  out  by  a  white-light  source2.  The 
advantages  of  the  proposed  white-light  signal 
processing  technique  are:  1.  It  is  capable  of 
suppressing  the  coherent  noise;  2,  White-light 
sources  ere  usually  Inexpensive;  3,  The  pro¬ 
cessing  environments  are  not  critical;  A.  The 


white-light  system  is  relatively  easy  and  econom¬ 
ical  to  maintain;  and  5.  The  white-light  pro¬ 
cessor  is  particularly  suitable  for  color  image 
processing. 

One  question  that  the  reader  may  ask,  since 
the  white-light  ayetem  offers  all  these  glamorous 
merits,  why  it  has  been  Ignored  for  so  long?  The 
answer  to  this  question  is  that,  it  was  a 
general  acceptance  that  an  incoherent  source 
cannot  process  the  signal  in  complex  amplitude. 
However,  none  of  the  practical  sources  are 
strictly  incoherent,  even  a  white-light  source. 

In  fact,  we  were  able  to  utilize  the  partial 
coherence  of  a  white-light  source  to  perform  the 
complex  amplitude  processing.  The  proposed 
white-light  processor,  on  one  hand  it  is  capable 
of  suppressing  the  coherent  noise  like  an  inco¬ 
herent  processor,  on  the  other  hand  it  is  capable 
of  processing  the  signal  in  complex  amplitude  like 
a  coherent  processor. 

There  is  however  a  basic  different  approach 
toward  a  coherent  and  a  white-light  procesaor. 

In  coherent  processing,  virtually  no  one  seems  to 
care  about  the  coherence  requirements,  since  the 
laser  provides  a  strong  coherent  source.  However, 
in  white-light  processing,  the  knowledge  of  the 
coherence  requirement  is  usually  needed. 

In  white-light  processing  we  would  approach 
the  problem  backward.  First,  we  should  know 
what  is  the  processing  operation  we  wish  to  per¬ 
form:  Is  it  a  1-D  or  2-0  processing?  Is  the 
signal  filtering  a  point  or  point-pair  concept? 
What  is  the  spatial  bandwidth  of  the  signal?  etc. 
Then  with  this  knowledge,  we  would  be  able  to 
evaluate  the  coherence  requirements  at  the  Fourier 
and  at  the  input  planes.  From  the  evaluated 
results,  we  would  be  able  to  design  a  signal 
sampling  function  and  a  source  encoding  function 
to  obtain  these  requirements.  The  objective  of 
using  a  signal  sampling  function  is  to  achieve  a 
high  dogroe  of  temporal  coherence  in  Fourier  plane 
so  that  the  signal  can  be  processing  in  complex 
amplitude,  for  the  entire  spectral  band  of  a 
white-light  source.  And  for  the  source  encoding 
is  to  alleviate  the  constraints  of  an  extended 
white-light  source. 

In  the  following  sections,  we  shall  discuss 
la  detail  the  source  encoding,  signal  sampling 
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and  {Hearing  aa  appliad  to  a  white-light  signal 
processing. 
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White  Light  Signal  Processing 

We  shall  now  describe  an  optical  signal  pro¬ 
cessing  technique  that  can  be  carried  out  by  a 
white-light  source,  as  Illustrated  in  Fig.  1.  The 
white-light  signal  processing  system  is  similar  to 
that  of  a  coherent  system,  except  the  use  of  a 
white-light  source,  source  encoding  mask,  signal 
sampling  grating,  and  achromatic  transform  lenses. 
For  example,  if  we  place  a  signal  transparency 
s(x,y)  in  contact  with  a  sampling  phase  grating, 
the  complex  light  field  for  every  wavelength  X 
behind  the  achromatic  transform  lens  would  be 

E(P,q;A)  “  //s(x,y)exp(ipox)expl-l(px+qy)]dxdy 

-  S(p-po,q),  (1) 

where  the  integral  is  over  the  spatial  domain  of 
the  input  plane  P.,  (p,q)  denotes  the  angular 
spatial  frequency  coordinate  system,  p  is  the 
angular  spatial  frequency  of  the  sampling  phase 
gracing,  and  S(p,q)  is  the  Fourier  spectrum  of 
s(x,y).  If  we  write  Eq.  (1)  in  the  form  of 
linear  spatial  coordinate  system  (a, 6),  we  have, 

E(a,6;X)  -  S(a  -  po,0),  (2) 

where  p  ^  (2ir/Xf)a,  q  ^  (2ir/Xf)B,  and  f  is  the 
focal  length  of  the  achromatic  transform  lens. 
Thus,  we  see  that  the  Fourier  spectra  would  dis¬ 
perse  into  rainbow  color  along  the  a  axis,  and 
each  Fourier  spectrum  for  a  given  wavelength  X 
is  centered  at  a  +  (Af/2ir)po. 


Fig.  1.  A  white-light  optical  signal  processor. 


In  signal  filtering,  we  assume  that  a 
aequence  of  complex  spatial  filters  for  various 

X.  are  available,  l.e.,  ll(p_,q_) ,  whore 
o  d  n 

Pn  -  (2«/Xnf)a,  qn  »  (2w/Xnf)6.  In  practice, 

all  the  processing  signals  are  spatial  frequency 
limited,  the  spatial  bandwidth  of  each  H(p  ,q) 
is  also  band 11ml tad,  such  as 


,  fH<W>  V  0<  °2* 

(_  0,  otherwise  » 


(3) 


where  -  (X^f /2rr)  <P0  +  A p)  and 

Oj  «  (Xnf/2xXpo  -  Ap)  are  the  upper  and  the 

lower  spatial  limits  of  H(p  , q  ),  and  Ap  is  the 

n  n 

spatial  bandwidth  of  the  input  signal  s(x,y). 


The  limiting  wavelengths  of  each 
be  written  as 

p  +  Ap  p  -  Ap 

A,  »  A  - t—  ,  and  A.  «  A  - r ' 

l  n  p  -  Ap  h  np+Ap 

o  o 

The  spectral  bandwidth  of  H(pQ,qn)  la 


H(p.q„)  can 
n  q 

(4) 


therefore. 


AX 

n 


X 


n 


4pQAp 

1  7k  n2 

P  -  (Ap) 


^ x  . 

po  " 


(5) 


If  we  place  this  set  of  spatial  filters  side-by- 
side  positioned  over  the  smeared  Fourier  spectra, 
then  the  intensity  distribution  of  the  output 
light  field  can  be  shown  as, 

N  ? 

I(x,y)  *  Z  AXn  |s(x,y;Xn)*h(x,y;Xn)r,  (6) 

n*l 


where  h(x,y;X)  is  the  spatial  impulse  response  of 
H(pa,qn)  and  *  denotes  the  convolution  operation. 

Thus,  the  proposed  white-light  signal  processor  Is 
capable  of  proceeslng  the  signal  in  complex  ampli¬ 
tude.  Since  the  output  intensity  Is  the  sum  of 
the  mutually  incoherent  narrow  band  lrradiancee, 
the  annoying  coherent  noise  can  be  eliminated. 
Furthermore,  the  white-light  source  contains  all 
the  color  wavelengths,  the  proposed  system  is 
particularly  suitable  for  color  signal  processing. 

Spatial  Filtering.  Signal  Sampling 
and  Source  Encoding 

We  have  mentioned  earlier  for  white-light 
processing,  we  would  approach  the  problem  in 
backward  manner.  For  example,  if  signal  filter¬ 
ing  is  two-dimensional  (e,g,,  2-0  corralation 
operation),  we  would  synthesize  a  set  of  narrow 
spectral  band  filters  for  each  XR  for  the  entire 
smeared  Fourier  spectra,  as  illustrated  in 
Fig,  2(a),  On  the  other  hand,  if  the  signal 
filtering  is  one-dimensional  (e.g.,  deblurring 
due  to  linear  motion) ,  a  broadband  fan-shape 
spatial  filter,  to  accomodate  the  scale  variation 
due  to  wavelength,  can  be  utilized  as  illustrated 
in  Fig.  2(b),  Since  the  filtering  is  taken  place 
with  the  entire  spectral  band  of  the  light  source, 
the  coherent  noise  can  be  suppressed  and  tha 
white-light  procasslng  technique  is  very  suitable 
for  color  Image  processing. 

There  is,  however,  a  temporal  coherence 
requirement  imposed  upon  the  signal  filtering  in 
Fourier  plane.  Since  the  scale  of  the  Fourier 
spectrum  varies  with  wavelength,  a  temporal 
coherence  requirement  should  be  Imposed  on  each 
■patlal  filter  at  tha  Fourier  plana.  Thus,  the 
spectral  spread  over  each  filter  H(p  ,q)  la 
Imposed  by  the  temporal  coherence  requirement , 
l.e., 
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(a)  A  multi-spectral-band  filter. 


0>)  A  fan-shape  filter 
Fig.  2. 


ME  «  1. 

Po 


(7) 


From  this  requirement ,  a  high  degree  of  temporal 
coherence  is  achievable  In  the  Fowler  plane  by 
simply  Increasing  the  spatial  frequency  of  the 
sampling  grating.  Heedless  to  say  that,  the 
same  temporal  coherence  requirement  of  Eq.  (7)  can 
also  be  applied  for  a  broadband  fan-shape  filter. 

There  la  also  a  spatial  coherence  require¬ 
ment  Imposed  at  the  Input  plane  of  the  whlte- 
llght  processor.  With  reference  to  the  Wolf's 
partial  chherence  theory-* ,  the  spatial  coherence 
function  at  the  input  plage  can  be  shown  . 

r(x-x')  -  //y(Jo)exp[i2»  ^  (Z-2')]<£q,  (8) 

where  y(x  )  denotes  the  Intensity  distribution  of 
o 

the  source  encoding  function. 

From  the  above  equation,  ve  see  that  the 
spatial  coherence  and  source  encoding  function 
form  a  Fourier  transform  pair,  l.e. , 

V(20)  -*£r<J-*')],  (9) 

and 

!•(*-*')  -  fl[Y<*#>],  (10) 

where  y  denotaa  the  Fourier  traasfermatlon.  This 

Fourier  transform  pair  Implies  that  If  a  spatial 
coherence  function  Is  given  then  tho  source 
encoding  function  can  be  evaluated  through  the 


Fourier  transformation  and  vice  versa.  We  note 
that  the  source  encoding  function  can  consist  of 
apertures  of  any  shape  or  complicated  gray  scale 
transparency.  However  the  source  encoding  func¬ 
tion  Is  only  limited  to  a  positive  real  quantity 
which  Is  restricted  by  the  following  physical 
realizable  condition: 

0iY(io)  il.  (11) 

In  white-light  processing,  we  would  search 
for  a  reduced  spatial  coherence  requirement  for 
the  processing  operation.  With  reference  to  this 
reduced  spatial  coherence  function,  a  source 
encoding  function  that  satisfied  the  physical 
realizability  condition  ca„  be  obtained.  One  of 
the  basic  objectives  of  the  source  encoding  Is  to 
alleviate  the  constraint  of  a  white-light  source. 
Furthermore  the  source  encoding  also  Improves  the 
utilization  of  the  light  source  such  that  the 
optical  processing  can  be  carried  out  by  an 
extended  source. 


We  shall  now  Illustrate  an  application  of 
the  source  encoding,  signal  sampling  and  filter¬ 
ing  for  a  white-light  signal  processing.  Let  . 
us  nov  consider  a  polychromatic  image  subtraction3. 
The  image  subtraction  of  Lee”  that  we  would  con¬ 
sider  Is  essentially  a  one-dlmenslonal  processing 
operation,  in  which  a  1-D  fan-shape  diffraction 
grating  should  be  utilized,  as  illustrated  In 
Fig.  3.  We  note  that  the  fan-shape  grating 
(l.e.,  filter)  Is  imposed  by  the  temporal  coher¬ 
ence  condition  of  Eq.  (7).  Since  the  Image  sub¬ 
traction  Is  a  point-pair  processing  operation, 
a  strictly  broad  spatial  coherence  function  at 
the  Input  plane  Is  not  required.  In  other  words, 
If  one  maintains  the  spatial  coherence  between  the 
corresponding  Image  points  to  be  subtracted  at 
the  input  plane,  then  the  subtraction  operation 
can  be  carried  out  at  the  output  Image  plane. 

Thus  instead  of  using  a  strictly  broad  spatial 
coherence  function,  a  reduced  spatial  cbbercnce 
function  may  be  utilized,  such  as 


r(y-y')  -  «(y-y'-ho)  +  «(y-y*+ho),  (12) 


where  2h  Is  the  main  separation  between  the  two 
Input  coSor  transparencies.  The  source  encoding 
function  can  therefore  be  evaluated  by  through 
the  Fourier  transform  of  Eq.  (9),  such  as 
2»h 

Y(y0)  -  2  cos  (-j^yo).  (13) 


Unfortunately  Eq.  (13)  Is  a  bipolar  function  which 
la  not  physically  realizable.  To  ensure  a  physi¬ 
cally  realizable  source  encoding  function,  we  let 
a  reduced  spatial  coherence  function  with  the 
required  point-pair  coherence  characteristic  be: 


r(lT-r’l) 


,ln<h1  ly"5r’l) 

_ o 

Hsln(|-  |y-y'|) 
o 


,lnc(FJd  Ijmt'I).  <W 

o 


where  H  »  1  a  positive  Integer,  and  w  «  d. 
Equation  (14)  represents  s  sequence  of  narrow 
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pulses  which  occur  aC  every  | y— y ’ |  •  nhQ,  where 

n  is  i  positive  integer,  and  their  peek  values 
are  weighted  by  a  broader  sine  factor,  as  shown 
in  Fig.  4(a).  Thus,  a  high  degree  of  spatial 
coherence  can  be  achieved  at  every  point-pair 
between  the  two  input  color  transparencies.  By 
taking  the  Fourier  tranof urination  of  the  reduced 
spatial  coherence  function  of  Eq.  (14),  the 
corresponding  source  encoding  function  is 

Y(|y|)  •  E  rect  >  (15) 

n-1  w 

where  w  is  the  slit  width,  d  -  (Xf/h  )  is  the 
separation  between  the  slits,  and  N  ?s  the  number 
of  the  slits.  Since  y(|y|)  is  a  positive  real 
function  which  satisfies  the  constraint  of 
Eq.  (11),  the  proposed  source  encoding  function 
of  Eq.  (15)  is  physically  realisable. 


Fig.  3,  A  white-light  image  subtraction 
processor.  TOO;  phase  grating, 

Lgl  image  Ians,  Lc;  collimated  lens, 

Lj  and  Lji  achromatic  transform 

lenses,  y(y)j  source  encoding  mask, 

G;  fan-shape  diffraction  grating. 

In  view  of  Eq.  (15),  we  also  note  that,  the 
separation  of  slit  d  is  linearly  proportional  of 
the  X.  The  source  encoding  is  a  fan-shape  type 
function,  as  shown  in  Fig,  4(b).  To  obtain  lines 
of  rainbow  color  spectral  light  source  for  the 
signal  processing,  we  would  utilize  a  linear 
extended  white-light  source  with  a  dispersive 
phase  grating,  as  illustrated  in  Fig,  3,  Thus 
with  appropriate  source  encoding,  signal  sampling 
and  filtering,  color  image  subtraction  operation 
can  be  obtained  at  the  output  plane.  We  stress 
again,  the  basic  advantage  of  source  encoding 
is  to  alleviate  the  constraint  of  strict  spatial 
coherence  requirement  imposed  upon  the  optical 
signal  processor.  The  sourca  encoding  also 
offers  the  advantage  of  efficient  utilisation 
of  the  light  power. 


n»-p'i 
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(b)  A  source  encoding  mask. 
Fig.  4. 


perlmental  Demonstrations 


We  shall  now  provide  a  couple  of  experimental 
results  obtained  with  the  source  encoding,  signal 
sampling  and  filtering  technique  for  white-light 
and  extended  incoherent  sources.  We  shall  first 
show  the  result  obtained  for  color  image  deblur¬ 
ring  due  to  linear  motion  with  the  white-light 
processing  technique.  Since  linear  motion 
deblurring  is  a  1-D  processing  operation  end 
the  inverse  filtering  is  a  point-by-point  filter¬ 
ing  concept  such  that  the  operation  is  taking 
place  on  the  smearing  length  of  the  blurred 
image.  Thus  the  deblurring  filter  (l.e.. 

Inverse  filter)  is  a  fan-shape  type  spatial 
filter”  and  the  temporal  coherence  requirement 
is  imposed  by  Eq,  (11).  The  spatial  coherence 
requirement  is  dependent  upon  the  smearing 
length,  A  source  encoding  function  of  a 
narrow  slit  width  (dependent  upon  the  smearing) 
perpendicular  to  the  smearing  length  is  utilized. 
Figure  5(e)  shows  a  black-and-white  color  picture 
of  s  blurred  image  due  to  linear  motion  of  a 
F-16  fighter  plane.  The  body  of  the  fighter 
plane  is  painted  in  black-and-white  colors,  the 
wings  are  mostly  painted  in  red,  the  tail  la 
black-and-white,  and  the  ground  terrain  is 
generally  bluish  color.  From  this  figure,  we 
see  that  the  plane  is  badly  blurred.  Figure  5(b) 
shows  the  color  image  deblurringr  result  that  we 
obtained  with  the  proposed  white-light  deblurring 
technique.  From  this  deblurrdd  result,  the 
letters  and  overall  shape  of  the  entire  airplane 
are  more  distinctive  than  the  blurred  one. 
Furthermore,  the  river,  the  highways,  and  the 
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forestry  of  the  ground  terrain  are  far  more 
visible.  We  note  that  the  color  reproduction  of 
the  deblurred  image  is  spectacularly  faithful, 
and  coherent  artifact  noise  is  virtually  non- 
existed.  There  is,  however,  some  degree  of  color 
blur  and  color  deviation,  which  are  primarily 
due  to  the  chromatic  aberration  and  the  anti- 
ref lcct.'inco  continp,  of  the  transform  lenses. 
Nevertheless,  these  drawbacks  can  be  overcome  by 
utilizing  good  quality  achromatic  transform 
lenses. 


(a)  A  black-and-white  picture  of  a  blurred 
color  Image. 


(b)  A  black-and-white  picture  of  a  de¬ 
blurred  color  image. 

Fig.  5. 

hot  us  now  provide  n  color  imago  subtrac¬ 
tion  utilized  by  the  source  encoding  technique 
with  extended  incoherent  sources  as  described  in 
previous  sections.  Figure  6(a)  and  6(b)  shows 
two  black-and-white  color  Image  transparencies 
of  a  parking  lot  as  Input  color  objects. 

Figure  6(c)  shove  the  color  subtracted  image 
obtained  by  the  source  encoding  technique  with 
extended  incoherent  source.  In  this  figure,  the 
profile  of  a  (red)  subcompact  car  can  be  seen  at 


the  output  image  plane.  The  shadow  and  the 
parking  line  (in  yellow  color)  can  also  be 
readily  identified.  We  however  note  that,  this 
color  image  subtraction  result  is  obtained  by  two 
narrow  band  extended  incoherent  sources.  Exten- 
tion  toward  the  entire  spectral  band  of  a  white- 
light  source  is  currently  under  investigation. 


(a)  A  black-and-white  picture  of 
the  input  color  objects. 


(b)  A  black-and-white  picture  of 
the  input  color  objects. 
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(c)  A  black-and-white  picture  of  the 
subtracted  color  image. 

Fig.  6. 

Conclusion 

In  conclusion  we  would  point  out  that  the 
advantage  of  source  encoding  is  to  provide  an 
appropriate  spatial  coherence  function  at  the 
input  plane  so  that  the  signal  processing  can 
be  carried  out  by  an  extended  incoherent  source. 
The  effect  of  the  signal  sampling  is  to  achieve 
the  temporal  coherence  requirement  at  the 
Fourier  plane  so  that  the  signal  can  be  processed 
in  complex  amplitude.  If  the  filtering  operation 
is  two-dimensional,  a  multi-spectral-band  2-0 
filters  should  be  utilized.  If  the  filtering 
-operation  is  one-diaiensional,  a  fan-shape  filter 
can  be  used. 

In  short,  one  should  carry  out  the  proces¬ 
sing  requirements  backward  for  a  white-light  pro¬ 
cessing.  With  these  processing  requirements 
(e.g.,  operation,  temporal  and  spatial  coherence 
requirements),  multi-spectral-band  or  broad-band 
filter,  signal  sampling  function,  and  source 
encoding  mask  can  be  synthesized.  Thus  the  signal 
processing  can  be  carried  out  in  complex  amplitude 
over  the  whole-spectral  band  of  the  light  source. 

We  acknowledge  the  support  of  the  0.S,  Air 
Force  Office  of  Scientific  Research  Grant 
AF0SR-81-0148, 
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Soustractioo  d’images  en  couleur 
utiHsant  des  sources  et endues  incoherentes 

SUMMARY  :  A  technique  of  color  image  subtraction  with  encoded 
extended  incoherent  sources  is  presented.  The  objective  of  the  source 
encoding  is  to  obtain  a  reduced  coherence  requirement  for  the  image 
subtraction  so  that  the  inherent  difficulty  of  obtaining  an  incoherent 
point  source  can  be  alleviated.  Experimental  results  for  the  color 
image  subtraction  obtained  with  the  incoherent  processing  techni¬ 
que  are  given.  Since  most  images  are  multi-colored,  this  color  image 
subtraction  technique  may  offer  a  broad  range  of  applications. 


R&SUM&  :  Nous  prtsentons  une  mtthode  de  soustraction  d'images 
en  couleur  utilisant  des  sources  incoherentes  codecs.  Le  but  du 
codage  est  d'obtenir  la  coherence  necessairc  pour  notre  experience 
avec  une  source  relativement  grande.  Nous  donnons  des  rEsulut s 
experimentaux.  La  plupart  des  images  etant  des  images  en  couleur, 
cette  technique  de  soustraction  peut  avoir  des  applications  diverses 


1.  —  INTRODUCTION 

One  of  the  most  interesting  segments  of  optical 
information  processing  must  be  the  image  subtrac¬ 
tion-  The  applications  may  be  of  value  in  earth  resource 
studies,  urban  development,  highway  planning,  land 
use,  surveillance,  automatic  tracking,  and  many 
others.  The  image  subtraction  may  also  be  applied 
to  video  communications  as  a  means  of  bandwidth 
compression.  For  example  it  is  only  necessary  to 
transmit  the  differences  between  the  two  images 
in  successive  cycles,  rather  than  the  entire  image  in 
each  cycle. 

There  are  several  techniques  available  for  image 
subtraction  which  can  be  found  in  a  review  paper  by 
Ebersole  [I].  However  most  of  the  optical  image 
subtraction  techniques  have  relied  on  a  coherent 
source  to  carry  out  the  subtraction  operation.  But  a 
coherent  optical  processing  system  is  plagued  with 
coherent  noise,  which  frequently  limits  its  processing 
capabilities. 


(*)  Vetting  scholar  from  Shanghai  Institute  of  Optics  and  Fine 
Mechanics,  Academia  Sinica  (China) 


We  have  in  a  recent  paper  [2]  demonstrated  a 
technique  of  image  subtraction  utilizing  an  extended 
incoherent  source.  Since  the  image  subtraction  is  a 
one-dimensional  processing  operation  and  the  spa¬ 
tial  coherence  requirement  for  the  subtraction  ope¬ 
ration  is  a  point-pair  coherence  requirement  it  is, 
therefore,  possible  to  encode  an  extended  source  to 
obtain  the  required  spatial  coherence.  In  a  previous 
paper  [3]  we  have  shown  a  Fourier  transform  rela¬ 
tionship  between  the  spatial  coherence  and  source 
intensity  distribution.  In  principle,  it  is  possible  to 
encode  an  extended  source  to  obtain  an  appropriate 
spatial  coherence  for  specific  information  processing 
operations. 

Strictly  speaking,  all  images  in  the  visible  wave¬ 
lengths,  which  includes  the  black-and-white  images, 
are  color  images.  Therefore  it  is  of  interest  for  us, 
in  this  paper,  to  extend  this  incoherent  processing 
technique  for  color  image  subtraction. 

n.  —  COLOR  IMAGE  SUBTRACTION 

We  will  now  describe  a  color  image  subtraction 
operation  with  an  encoded  extended  incoherent 
source,  as  depicted  in  figure  I.  For  simplicity,  two 
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incoherent  light  sources,  each  for  a  different  color  of 
light  (<.«.,  red  and  green),  are  used  for  the  subtraction 
operation.  For  the  detailed  analysis  of  image  sub¬ 
traction  with  an  encoded  extended  source  we  refer 
the  reader  to  our  previous  paper  [2]. 

In  color  image  subtraction,  we  insert  the  color 
image  transparencies  Ov(x,  y)  and  02(jc,  y)  in  the 
open  apertures  of  the  input  plane  P2,  which  can  be 
described  as 

(1)  /(*.  y)  =  Oi(*  -  h0,  y)  +  02(x  +  h0,  y) , 

where  2  h0  is  the  separation  between  the  two  input 
transparencies  O,  and  02.  Two  sinusoidal  gratings 
G,  and  G2  designed  for  the  red  and  the  green  color 
wavelengths  respectively  are  inserted  in  the  spatial 
frequency  plane  P„  and  can  be  written  as 

(2)  G,  «  ^  [1  +  sin  (A0Pr)3  ■ 
and 

(3)  G2  =i[l  +  an  (ft,*)], 

where  p,  «  (2  *«)/(*,  f)  and  pt  -  (2  na)HXgf)  are 
the  spatial  frequencies  of  the  gratings,  X,  and  X,  are 
the  red  and  green  color  wavelengths,  a  denotes  the 
spatial  coordinate  in  the  same  direction  as  p,  and  /  is 
the  focal  length  of  the  achromatic  transform  lens  L} 
By  a  straight-forward  but  rather  cumbersome  eva¬ 
luation,  the  irradiance  around  the  origin  of  the  output 
image  plane  Pt  can  be  shown  [2], 

(4)  A*.  y)  -  /,(*,  y)  +  ijix,  y) 

-  K  |  0„(x,  y)  -  02r(x,  y)  |J  + 

+  K  |  0,/x,  y)  -  02jx,  y )  |* , 


where  /,(*,  y)  and  Ij(x,  y)  denote  the  red  and  green 
subtracted  color  image  irradiances,  Olr,  Oit 
and  02f  are  the  corresponding  red  and  green  color 
input  oojects.  From  Eq.  (4)  we  see  that  the  subtract¬ 
ed  color  image  can  be  obtained  at  the  output  image 
plane-  Since  the  color  image  subtraction  is  obtained 
with  extended  incoherent  sources,  the  coherent  arti¬ 
fact  noise  can  be  suppressed. 

in.  —  SOURCE  ENCODING 
FOR  IMAGE  SUBTRACTION 

We  would  now  consider  a  source  encoding  technique 
for  image  subtraction.  As  we  have  pointed  out  earlier, 
the  image  subtraction  is  a  1  -  D  processing  opera¬ 
tion  and  the  spatial  coherence  requirement  is  depen¬ 
dent  upon  corresponding  image  points  to  be  subtract¬ 
ed.  Therefore,  instead  of  utilizing  a  strictly  broad  spa¬ 
tial  coherence  function,  a  point-pair  spatial  coherence 
is  sufficient.  To  insure  a  physically  realizable  source 
intensity  distribution,  we  let  the  point-pair  spatial 
coherence  be  [2,  3] 


where  2  h0  is  the  separation  between  the  image  points, 
N  ►  1  a  positive  integer,  and  w  <  d  The  first  factor 
of  Eq.  (3)  represents  a  sequence  of  narrow  pulses 
present  at  every  |  x  -  x'  \  nh *  where  n  is  a  positive 
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integer,  and  the  peak  values  of  the  pulses  are  weighted 
by  a  broader  sine  factor,  as  shown  in  figure  2(a). 
From  this  figure,  we  see  that  there  exists  a  high  degree 
of  spatial  coherence  between  every  point-pair  of  the 
two  input  objects  Ox  and  02-  By  taking  the  Fourier 
transformation  of  Eq.  (5),  we  obtain  the  following 
intensity  distribution  [3] 

(6)  5(|  x  |)  =  £  rect  ( '*—■—!  , 

where  w  is  the  slit  width,  d  =  (A f)/h0  is  the  separation 
between  the  slits  of  the  source  coding  apertures,  A 
is  the  wavelength  of  the  light  source,  and  /  is  the  focal 
length  of  the  collimating  lens.  Thus  the  intensity 
distribution  of  Eq.  (6)  represents  N  number  of  narrow 
slit  apertures  with  equal  spacing  d,  as  shown  in 
figure  2(b).  From  this  result,  it  is  possible  to  encode 
an  extended  source  to  obtain  a  spatial  coherence  func¬ 
tion  at  the  input  plane  for  the  subtraction  operation. 


HI - If— 
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Fig.  2.  —  (»)  A  Point-pair  Spatial  Coherence  Function,  (b)  Source 
Encoding  Mask. 


Since  the  scale  of  the  Fourier  spectrum  varies  with 
the  wavelength  of  the  light  source,  a  temporal  cohe¬ 
rence  requirement  should  also  be  imposed  on  every 
processing  operation.  In  other  words,  the  spectral 
width  (i.e.,  temporal  coherence)  of  the  light  source 
should  be  restricted  by  the  following  inequality  [3] 


(7) 


M  _n_ 
*■  *  *0  Pm 


spatial  frequency  requirement  is  high,  then  a  narrow 
spectral  width  of  the  light  source  is  required. 

IV.  —  EXPERIMENTAL  DEMONSTRATION 

Strictly  speaking,  light  sources  emitting  all  primary 
colors  (i.e.,  red,  green,  and  blue)  should  be  used  for 
the  color  images  subtraction.  For  simplicity  of  expe¬ 
rimental  demonstration,  a  mercury  arc  lamp  with  a 
green  filter  (5  461  A)  and  a  zirconium  arc  lamp  with 
a  red  filter  (6  328  A)  were  used  for  the  color  light 
sources.  The  intensity  ratio  of  the  two  resulting  light 
sources  was  adjusted  to  about  unity  with  a  variable 
beam  splitter. 

The  slit  widths  for  the  source  encoding  masks  were 
about  2.5  p  and  the  spacings  of  the  slits  were  25  p 
for  the  green  wavelength  and  29  p  for  the  red  wave¬ 
length.  The  overall  size  of  the  source  encoding  masks 
was  about  3x3  mm2  which  contained  about  120 
and  100  slits  respectively.  The  focal  length  of  the 
transform  lenses  was  300  mm.  A  liquid  gate  containing 
two  color  image  transparencies  about  6x8  mm2 
each  and  with  a  separation  of  about  13.2  mm  was 
placed  behind  the  collimator.  Two  sinusoidal  gratings 
with  spatial  frequencies  of  1/(25  p)  and  1/(29  p)  were 
used  for  the  green  and  red  color  image  subtraction 
operation,  as  shown  in  figure  I. 

In  our  first  experimental  demonstration,  we  provide 
in  figures  3(a)  and  3(b)  two  sets  of  different  colored 


where  A  is  the  center  wavelength  of  the  light  source 
and  pm  is  the  highest  angular  spatial  frequency  requir¬ 
ed  for  the  image  subtraction.  Consequently,  if  the 


Fio.  3.  —  Color  Image  Subtraction ;  Binary  Object,  (a)  et  (b)  Input 
Color  English  Words,  (c)  Subtracted  Color  Image 
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English  words  as  input  objects.  Figure  3(a  shows 
the  subtracted  color  image  obtained  with  the  color 
subtraction  technique.  The  input  words  STATE 
L'S’IV.  are  red  and  green  respectively  in  one  set  but 
they  are  green  and  red  respectively  in  the  other  set 
while  in  the  subtracted  image,  as  shown  in  figure  3  (a, 
both  words  are  in  yellow  color,  which  is  consistent 
with  the  result  we  expected.  In  other  words,  the  sub¬ 
tracted  image  of  red  and  green  produces  yellow  color 
since  the  red  and  green  wavelengths  are  incoherent 
and  therefore  add  incoherently  to  produce  a  yellow 
color. 

For  a  second  demonstration,  we  provide  two  conti¬ 
nuous  tone  color  images  of  two  sets  of  fruit,  as  shown 
in  figures  4(a)  and  4(b).  By  comparing  these  two 
figures,  we  see  that  a  dark  green  cucumber  and  a  red 
tomato  are  missing  in  figure  4(b).  Figure  4(c)  shows 
the  subtracted  color  image  obtained  with  this  incohe¬ 
rent  color  image  subtraction  technique.  In  this  result. 
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the  profiles  of  a  cucumber  and  a  tomato  can  be  seen 
at  the  output  image  plane. 

For  the  final  demonstration,  we  again  provide 
two  continuous  tone  color  images  of  a  parking  lot  as 
input  color  object  transparencies,  as  shown  in 
figures  5(a)  and  5(b).  From  these  input  transparen¬ 
cies.  we  see  that  a  red  color  passenger  car  shown  in 
the  parking  lot  in  figure  5(a)  is  missing  in  figure  5(bj. 
Figure  5(c)  is  the  color  subtracted  image  obtained 
from  the  incoherent  color  image  subtraction  techni¬ 
que  as  previously  described.  In  this  figure,  a  red 
passenger  car  can  clearly  be  seen  at  the  output  image 
plane,  lt  is  also  interesting  to  point  out  that  the  parking 
line  (in  yellow  color)  on  the  right  side  of  the  red  car 
can  readily  be  seen  with  the  subtracted  image. 


Fio.  4.  —  Color  Image  Subtraction .  Continuous  Tone  Object 
(a)  et  (b)  Input  Color  Object  Transparencies  Subtracted  Color 
Image. 


Flo.  5  Color  Image  Subtraction .  Parking  Lot  tal  et  (bl  Input 
Color  Transparencies  Icl  Subtracted  Color  Image 
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V.  —  CONCLUSION 

We  have  introduced  a  color  image  subtraction 
technique  with  encoded  incoherent  sources.  The 
bask  advantage  of  source  encoding  is  to  increase 
the  available  light  power  for  the  image  subtraction 
operation,  so  that  the  inherent  difficulty  of  obtaining 
incoherent  point  sources  can  be  alleviated.  Since  the 
technique  uses  incoherent  sources,  the  annoying 
coherent  artifact  noise  can  be  suppressed.  We  would 
see  that  the  concept  of  color  image  subtraction  may 
also  be  extended  to  the  use  of  white-light  source,  for 
which  a  program  is  currently  under  investigation.  In 
experimental  demonstrations,  we  have  shown  that 
color  subtracted  images  can  be  easily  obtained  by  this 


incoherent  subtraction  technique.  Since  virtually  all 
ia sage*  of  natural  objects  are  color,  the  technique  may 
offer  a  wide  range  of  practical  applkations. 
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Rainbow  holographic  aberrations  and  the  bandwidth 
requirements 

Y.  W.  Zhang,  W.  G.  Zhu.  axl  F.  T.  S.  Yu 

Rainbow  holographic  image  resolution,  primary  aberrations,  and  bandwidth  requirements  are  presented. 
The  results  obtained  for  the  rainbow  holographic  process  are  rather  general,  for  which  the  conventional  holo¬ 
graphic  image  resolution,  aberrations,  and  bandwidth  requirements,  can  be  derived.  The  conditions  for  the 
elimination  of  the  five  primary  rainbow  holographic  aberrations  are  also  given.  These  conditions  may  be 
useful  for  the  application  of  obtaining  a  high-quality  rainbow  hologram  image.  In  terms  of  bandwidth  re¬ 
quirements,  we  have  shown  that  the  bandwidth  requirement  for  a  rainbow  holographic  construction  it  usually 
several  orders  lower  than  that  of  a  conventional  holographic  process.  Therefore,  a  lower-resolution  record¬ 
ing  medium  can  generally  be  used  for  most  of  the  rainbow  holographic  constructions. 


I.  Introduction 

The  rainbow  holographic  process  of  Benton1  involves 
two  recording  steps.  First,  a  primary  hologram  is  made 
from  a  real  object  with  the  conventional  off-axis  holo¬ 
graphic  techniques,  and  second,  the  rainbow  hologram 
is  recorded  from  the  real  hologram  image  from  the  pri¬ 
mary  hologram.  The  relaxation  of  the  coherence  re¬ 
quirement  on  the  reconstructing  process  arises  from  the 
placement  of  a  narrow  slit  behind  the  primary  hologram 
in  the  second-step  holographic  recording.  However,  the 
two-step  holographic  recording  process  is  cumbersome 
and  requires  a  separate  optical  setup  for  each  step.  It 
is  usually  a  major  undertaking  for  laboratories  with 
limited  optical  components.  An  alternative  method  of 
obtaining  a  rainbow  hologram  is  the  one-step  process 
of  Chen  et  al.2~ 5  They  had  shown  that  the  one-step 
technique  offers  certain  flexibilities,  and  the  optical 
arrangement  is  simpler  than  the  conventional  two-step 
process.  We  also  note  that  the  color  blur  of  the  Benton 
type  rainbow  hologram  has  been  subsequently  analyzed 
by  Wyant,6  Chen,7  and  Tamura.8  Recently,  Zhuang  et 
al. 9  investigated  the  image  resolution  and  color  blur  for 
a  one-step  rainbow  holographic  process  utilizing  dif¬ 
fraction  optics. 

In  this  paper,  we  shall  evaluate  the  primary  aberra¬ 
tions  and  bandwidth  requirement  for  a  rainbow  holo¬ 
graphic  process.  The  correlations  for  the  elimination 
of  the  rainbow  holographic  aberrations  are  presented. 
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Since  the  effective  aperture  of  a  rainbow  hologram 
is  generally  smaller  than  that  of  a  conventional  holo¬ 
gram,  the  spatial  bandwidth  requirement  is  generally 
somewhat  lower  for  the  rainbow  holographic  process. 
Therefore,  it  is  usually  possible  to  construct  a  rainbow 
hologram  utilizing  a  lower-resolution  Him. 

RMwOw  novogrspnic  hotovuhoh 

In  a  recent  paper,6  the  resolution  limit  and  color  blur 
of  a  rainbow  hologram  were  analyzed  from  diffraction 
optics.  However,  in  practice,  the  slit  used  in  a  rainbow 
holographic  process  is  many  orders  larger  than  the  light 
wavelength.  Therefore,  the  rainbow  holographic  pro¬ 
cess  can  also  be  approached  from  slightly  different 
configurations  as  depicted  in  Fig.  1.  For  simplicity, 
both  the  recording  and  the  reconstruction  processes  are 
drawn  in  a  composite  diagram,  where  H \  represents  the 
primary  hologram  which  is  in  contact  with  a  narrow-slit 
aperture  SL,  H 2  represents  the  constructing  rainbow 
hologram,  O(xoJo) is the  object  image  point,  /(x,  o'i) » 
the  rainbow  hologram  image  point,  and  R(xryr)  and 
C(xciyc)  are  the  reference  and  the  reconstruction  point 
sources,  respectively.  We  note  that  if  O  and  SL  rep¬ 
resent  the  object  and  slit  images  due  to  an  imaging  lens, 
the  same  diagram  of  Fig.  1  can  also  be  used  for  the 
analysis  of  a  one-step  rainbow  holographic  process. 

With  reference  to  Fig.  1,  it  is  evident  that  the  coor¬ 
dinate  points  P „  and  P12  are  the  marginal  extensions 
of  the  rainbow  holographic  recording  region,  which  can 
be  written  as 


where  g  A  s  +  d. 
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Thus,  the  complex  light-amplitude  distribution  due 
to  the  object  image  point  O  at  the  rainbow  holographic 
recording  plate  H  2  is 

Ui(x\k\)  -  Ai  exp|-i*i[d2  +  (x  -  xo)2|1/2l,  (3) 

and  the  corresponding  complex  light  distribution  due 
to  the  reference  and  reconstructing  (called  illuminating) 
beams  are 


u2(x  jt ,)  -  AsexpHJkdi?  +  (X  -  x,)2l>'2|,  (4) 

U3(x-M)  -  A3  expt-i*2(X.|  +  (x  -  xc)2]>'2|,  (5) 

respectively,  where  the  spatial  impulse  response  is 

hi(a  — xMi)  *  At  expji*2(/2  +  (a  -  x)2]1/2|,  (6) 

Ai  is  a  complex  constant,  k  »  2r/\,  and  Ai  and  X2  are 
the  wavelengths  for  the  rainbow  hologram  construction 
and  the  reconstructing  processing,  respectively. 

The  corresponding  complex  light  field  of  the  rainbow 


hologram  image  point  I  is 

BiUrM  -  Ci 

« 

j*P\\  f  dw  1 

|  UiUtU^iidx  “  C2  sine  - — -  (a  +  a;)  • 

JPn  lAtft  1 

(7) 

where 

j  dLiLi 

ji/rili2  +  dL  1  —  fidL  2 

(8) 

/  *0  JCr  Xr\ 

(9) 

(/;,<*,)  identifies  the  position  of  the  rainbow  hologram 
image  point  /,  and  n  A  A2A1. 

The  magnification  of  the  hologram  image  can  then 
be  written  as 


Mi 


x0  d 


hZ+jr)" 


(10) 


and  the  resolution  limit  of  the  rainbow  holographic 
process  is 


AH, 


_ iLjid  ♦  s ) 

(fiL\L2  +  l>\d  “  itLid)W 


(ID 


where  IV  is  the  slit  width.  If  both  the  reference  and  the 
reconstruction  beams  are  plane  waves,  the  resolution 
limit  becomes 


f  « 


Fig.  1.  A  composite  rainbow  holographic  construction  and  recon¬ 
struction  process  for  the  evaluation  of  the  hologram  image  resolution 
and  the  primary  aberrations:  Hi,  primary  hologram;  SL,  slit  aper¬ 
ture;  W,  slit  width;  H2,  rainbow  holographic  plate;  O,  object  image 
point;  /,  rainbow  hologram  image  point;  ft,  convergent  reference  point 
source;  C,  white-light  reconstruction  point  source. 


AH, 


\,(d  +  s) 

W 


(12) 


We  note  that  these  results  are  identical  to  the  results  we 
obtained  previously  by  Zhuang  et  al. 10  It  is  evident 
that  this  partial  geometries  approach  simplifies  the 
rainbow  holographic  analysis. 


Hi.  Rainbow  Holographic  Aberrations 

In  evaluating  the  rainbow  holographic  aberration,  we 
expand  the  exponent  of  Eq.  (7)  in  binomial  expansion 
such  as 

v>(x»  -  *[£2  +  (x  -  o)2!,/2  -  *  |i  + <J  2L°>  -  +  J  • 

(13) 

Thus,  by  retaining  the  first  three  terms  of  Eq.  (13),  the 
third-order  primary  aberrations  of  a  rainbow  holo¬ 
graphic  process  can  be  obtained. 

To  evaluate  the  third-order  primary  aberrations,  the 
phase  factor  of  the  real  rainbow  hologram  image  can  be 
written  as 


Ifii  ■  *2Im(L?  +  <x  -  x,)2  +  (y  -  yr)2ll/2  -  p|d2  +  (x  -  x0)2  +  (y  -  y0)2],/2 

-  \L\  +  (x  -  x,)2  +  <y  -  yc)2!1'2  +  [/?  +  (a,  -x)2  +  (ft  - y)*|««|.  (14) 

With  reference  to  the  binomial  expansion  of  Eq.  (13), 
Eq.  (14)  can  be  written  as  a  configuration  of  the  order 
terms,  i.e., 


<ti  "  +  ft  +  >fh,  (15) 

where  <pc,  <t>P,  and  <t>h  are  the  first-,  second-,  and 
higher-order  terms,  which  shall  be  ignored  in  our  eval¬ 
uation.  The  third-order  term,  ,  is  written  as 

*>  -  j  l<«  -  *o)2  -My  -  yo)2|2  -  |U  -  xr)2 

+  (y  -  yr  >2]2  +  j-,  |(x  -  xc  >2  +  (y  -  y,  )2|2  -  ^  |(or,  -  x  )2  +  (ft  -  y )2)«j  •  (16) 
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Name 

Formula 

Condition  of  eliminating  aberration 

Spherical 

M  J*  .  1.  1 

I\  ds  If  L\ 

jj  “  l;<f  “  LiLi  m  L% 

c,  L,X r  Li 

ft «  l;d  •  U\L\  ~Lt~  «, 

Coma 

r  _  M  „  M  .  ft 
c>  L,yr  d3yo  +  t,  Li 

but,^ 1  7^  •  —■  •  jL  are  Finite 
hi  Li  L  2  ^2 

Astigmatism 

jut?  jut  J  o?  *? 

‘  Li  ~  d*  +  It  ~  Li 
.  jut  ry,  jutoyo  a, ft  Jtcy, 

'y  £?  ”  d2  .  /?,  “  Li 
.  .My,  j*y?.  0?  y? 
y  L?  “  d2  +  1?  ‘  L| 

jj  -  l;d  •  lit  ti  »  Lj  «  ®, 

but,  7^  >  7^  •  7s-  •  ^  are  finite 

4*1  4*i  4/2  4*2 

Curvature  of  field 

f-^(tt?  +  y?>-^(xl  +  yl) 

jj  »  l;  d  ■  1<;  Li  *  Lj  ■  «, 

+  i(a?  +  |3?)-^(x2+y2) 

but,  T-'T-’T-’T-  finite 

Distortion 

'  L,xr  d3*8+J?  t| 

n  _JLvs_Jlv}  +  2i_2£ 
y  L\y'  d2*S  i?  Li 

M  L,  Lj’L,  L* 

By  substituting  the  relation  of  x  ■  r  coed  and  y  ■  r  sinfl 
into  Eq.  (16),  die  third-order  phase  factor  can  be  shown 
as 


+  4  [d^  ^  +  y*>  “  fj (X'  +  y*>  +  Li <X*  +  y‘ > 


-  *,y«-  +  *«y«  -  fi  ajftJ  sin#  corfj  r2 

‘Lix’+Zixi-r?af}eo"+fcyi 


i 

2llw» 


-n>?  +  r;y?-^?|  •to#!' 


t!"  tr*  if 


(17) 


Thus,  can  be  written  as  a  function  of  the  five  primary 
aberrations10  such  as 


(18) 


where  S  represents  the  spherical  aberrations,  Cx  and 
Cy  represent  the  meridian  and  tangential  coma,  Ax,  Axy, 
and  Ay  are  the  astigmatisms,  F  is  die  curvature  of  field, 
and  Dm  and  Dy  are  the  horizontal  and  vertical  distor¬ 
tions.  By  comparing  this  equation  with  Eq.  (17),  the 
equations  for  the  five  primary  rainbow  holographic 
aberrations  can  be  obtained  as  tabulated  in  Table  I. 

In  this  table,  the  conditions  for  the  elimination  of  the 
primary  aberrations  are  given.  It  is  evident  that  all  the 
primary  aberrations  can  be  eliminated  if  both  the  ref¬ 
erence  and  the  reconstruction  beams  are  plane  waves 
andp«I.  It  is  interesting  to  note  that,  if  the  slit  width 
W  is  equal  to  the  size  of  the  primary  hologram,  by  let¬ 
ting  S  »  0,  the  primary  aberrations  reduce  to  the  con¬ 
ventional  off-axis  holographic  aberrations  as  obtained 
by  Meier.11 

Since  the  rainbow  holographic  imaging  is  generally 
produced  by  a  white-light  source,  the  aberrations  are 


ifit  “  *2  |~  ~  (C,  coo#  +  Cy  sin#)  —  —  F 


- (A,  co*2#  +  2 A„  sin#  co*#  +  A,  sin2#) 

2 


+  -  (D,  co*#  +  Dy  sin#) 
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affected  by  the  wavelength  spread  AA2  during  the  re¬ 
construction  process.  We  note  that  this  wavelength 
spread  AAS  is  limited  by  observer  eye  through  the  slit 
image.  Thus,  the  effect  of  aberration  variation  due  to 
AAj  can  be  determined  as  shown  in  the  following,  for 
simplicity  as  a  1-D  notation: 


From  Eq.  (21)  we  see  that  the  variation  of  the  astig¬ 
matism  is  identical  to  that  of  the  variation  of  curvature 
of  Held.  In  ending  this  section,  we  point  out  that  the 
above  equations  have  useful  applications  for  evaluating 
the  rainbow  hologram  image  quality. 


--£■*£& -rj* 

\r  mh  4A,  f  1  /*0  xrl 

‘  I1-  x.  ~  U* _  zfj  xT  Li?  It  ~  zj 


tjjti  —  jlpiLiLi  +  Lid  —  iiLyd)* 
(dLiLjt* 
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*0 
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Li  *L„ 
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_I| 

i/f£  _  ii) 

A,  C‘  +  A,  1 

m  dV 

ifl 

Id  LJ 

2 LtiLt-d) 


iltLiLt  +  Lid  —  nLfd ) 


(*o  .  *'  X, 
ltd  Li~>lLil 


"  •  “•  •(“  ir)lF  *  (fi  ■  si  tt  *  h1'  "n  ('?  *  8  '*fj*  !W-‘’  * w  (« ? 
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-  M7  +  f 

d  L  2 


dLiLt 


(19) 


(20) 


(21) 


(22) 


Fit-  X  A  composite  rainbow  boiofraphic  construction  and  racon- 
•traction  pnwaas  Ibr  dstennining  the  bandwidth  roquimnanta:  Hi, 
priaury  hokfram;  SL,  slit  apsrtun;  IV,  slit  width;  Hi,  rainbow  ho¬ 
lographic  plate;  Oi-Oj,  axtendad  hologram  object  imaca;  ft,  con  var¬ 
iant  rafaranea  point  aonroc. 


lea  BflnQWKRII  nv^UWliWiVlw 

Since  the  effective  aperture  of  a  rainbow  hologram 
is  generally  smaller  than  that  of  a  conventional  holo¬ 
gram,  it  is  possible  that  a  lower-resolution  film  can  be 
utilized.  Therefore,  it  is  our  aim  in  this  section  to  dis¬ 
cuss  the  bandwidth  requirements  of  a  rainbow  holo¬ 
graphic  process. 

We  shall  investigate  for  the  case  in  which  the  slit  and 
the  object  image  are  located  at  opposite  sides  of  the 
rainbow  holographic  plots  /f2  as  shown  in  Fig.  2,  where 
0 1  and  Oj  are  the  two  external  points  of  an  extended 
object  image  reconstructed  from  a  primary  hologram 
Hi.  Thus,  from  Fig.  2,  we  see  that 


Mi-3— iw- 

P,,-(l-^x„  +  ^IV. 

Pa“(i~3*°*"£w' 


(23) 
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Similar  to  what  we  did  in  previous  sections,  the 
complex  light  fields  U i,  I/2,  and  I/3  due  to  the  object 
image  points  0t  and  O2  and  the  reference  beam  R  can 
be  formulated.  The  terms  contributing  to  the  rainbow 
holographic  reconstruction  process  are 

(7,1/;  +  U3U[  -  2M,||A3|  CO.  {*,  |l,  -  d  -  -  ~^0l)2 

♦ttH- 

IhV 5  +  U3 V\  -  2|A2||A3|  CO.  {*,  ji.-d- 

where  A  2,  and  ^3  are  the  complex  constants,  and  8\ 
and  02  are  the  constant  phase  factors.  The  phase  shifts 
due  to  the  t/t,  U2,  and  I/3  are 

.  I.  .  U-*Ol)2.  U-XrW  .  „ 

- —  +  -^7“ )  +  <»•.  <»> 

+  +  (27) 

1  2d  2Li  ] 

respectively. 

The  derivatives  of  <t>3 1  and  <t> 32  represent  the  angular 
spatial  frequency  along  the  x  direction.  Thus,  the 
corresponding  linear  spatial  frequencies  are 


_ d*>31  m  J_  (M  1\  [  *01  srl 

2*  dx  XilUi  d)X  d  Lil 

_ d<P32  m  J_ ||J_ +  _ 

2x  dx  Xi  lUi  dl  d  Li] 


(28) 

(29) 


where  v  denotes  the  spatial  frequency.  Therefore,  the 
spatial  frequency  bandwidth  of  the  rainbow  hologram 
can  be  shown  as 


A»  "  *’32|j-/>n  -  m •  — —  |(£i  +  t)A xq 
Xi  Lig 


+  (L,-d)Wj. 

(30) 

where  g  *  s  +  d  and  Ajc0  * 

X02  ~  *oi-  If  we  let 

A*.  »  —  (l  +  — )  Axn, 

(31) 

\lg\ 

Lil 

A»2  "  r— |l 

-fW. 

(32) 

Xrf  \ 

Lil 

Eq.  (30)  reduces  to 

Ax  •  Axt 

+  Axit 

(33) 

where  Avi  and  Ai»2  are  the  spatial  frequency  band  widths 
due  to  the  size  of  object  image  and  the  slit  of  the  rainbow 
holographic  construction,  respectively.  Moreover,  Eq. 
(30)  can  also  be  written  as 

A*  *  ^  |ax0  +  VV  +  (sAx0-  dW)^-J  '  (34) 

From  this  equation,  we  see  that  if  (Ax0)/W  >  d/s,  in¬ 
creasing  L\  can  compress  the  spatial  frequency  band¬ 
width  of  the  rainbow  hologram.  On  the  other  hand,  if 
(Axo)/W  <  d/s,  decreasing  L\  can  also  compress  the 
spatial  frequency  bandwidth  of  the  rainbow  hologram. 
However,  for  the  case  (Ax0)/W  «  d/s,  the  L 1  is  inde¬ 
pendent  of  Av. 

Table  II  shows  the  spatial  frequency  bandwidth  re¬ 
quirements  for  various  rainbow  holographic  recording 
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conditions.  In  this  table,  the  notation  p  denotes  the 
distance  from  rainbow  holographic  plate  H 2  with  re¬ 
spect  to  the  reference  beam  R,  where  the  paraxial  ap¬ 
proximation  holds. 

We  note  that,  for  a  2-D  case,  the  general  formulas  for 
spatial  frequency  bandwidth  requirements  for  a  rain¬ 
bow  holographic  process  can  be  written  as 


1  1 

!  s 

d  1 

Ar,  -  - — 

Axo  + 

1  -  —  W 

1  L  i 

4il 

1  i 

1  8 

Ayo  + 

d\  1 

A",  «  - — 

1  +7- 

1-7-4 

1  Li 

4i|  J 

where  Avx  and  Avy  are  the  corresponding  spatial  fre¬ 
quency  bandwidths  in  the  x  and  y  directions,  respec¬ 
tively,  for  AXo  >  O,  Ayo  >  0,  and  s  >  0.  It  is  inter¬ 
esting  to  note  that  if  s  -  0,W  -  Lx,  and  L  -  Ly,  which 
correspond  to  the  recording  conditions  of  a  conventional 
hologram,  Eq.  (35)  reduces  to  the  results  obtained  by 
Yu,12  i.e., 


A*, 

Ar, 


l 
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X|d 


Axo  + 

i-ri 

l 

|Ay0  + 

>-fl 

I'll 

(36) 


To  have  some  feeling  of  magnitude,  we  would  provide 
two  numerical  examples  as  given  in  the  following: 

(1)  For  conventional  holograph  construction  process; 
we  let  Lx  *  Ly  *  100  mm,  Ax0  “  Ayo  ”  60  mm,  d  *  200 
mm,  -*  <»,  and  \i  ■  5  X  10-4  mm.  By  substituting 
these  data  into  Eq.  (36),  we  obtain  the  following  band¬ 
width  requirements:  Av'x  ■  Avy  “  1600  lines/mm. 

(2)  For  rainbow  holographic  construction  process, 
we  let  W  =  7  mm,  L  *  100  mm,  s  «  830  mm,  d  *  40  mm, 
L\  -  1000  mm,  Xx  *  5  X  10-4  mm,  and  Ax0  *  Ayo  *  6 
mm.  By  substituting  these  data  into  Eq.  (35),  we  obtain 
the  following  rainbow  holographic  bandwidth  require¬ 
ments  for  the  x  and  y  directions:  Avx  *  40  lines/mm, 
and  Auy  -  246  lines/mm.  In  comparison  with  these  two 
numerical  examples,  we  see  that  the  bandwidth  re¬ 
quirements  for  the  rainbow  holographic  construction 
are  somewhat  lower  than  those  of  the  conventional 
holographic  process.  Thus,  in  practice,  a  lower-reso¬ 
lution  film  can  be  utilized  for  rainbow  hologram  re¬ 
cording. 


V.  Conclusion 

The  image  resolution  and  third -order  primary  aber¬ 
rations  of  a  rainbow  holographic  process  were  derived. 
The  results  that  we  have  obtained  are  rather  general, 
for  which  the  resolution  limits  and  the  aberrations  of 
a  conventional  holographic  process  can  be  derived  from 
these  results.  The  conditions  for  the  elimination  of  the 
rainbow  holographic  aberrations  were  also  presented. 
These  conditions  are  very  useful  for  controlling  the 
rainbow  holographic  image  quality. 

The  bandwidth  requirement  of  a  rainbow  holographic 
process  was  also  evaluated.  As  compared  with  the 
conventional  holographic  process,  the  bandwidth  re¬ 
quirement  for  a  rainbow  holographic  process  is  usually 
several  orders  lower  than  that  of  a  conventional  holo¬ 
graphic  process.  Thus,  in  practice,  a  lower-resolution 
film  can  be  utilized  for  most  rainbow  holographic  con¬ 
structions. 
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